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ABSTRACT 


Research  reported  herein  is  directed  toward  further  understanding  of  the  fundamental 
processes  in  supersonic  combustion.  Part  One  describes  a  branching  chain  mechanism 
constructed  for  oxy-diborane  mixtures  diluted  in  argon.  Included  in  this  postulated 
mechanism  is  the  production  of  hydroxyl  as  an  ignition  intermediate.  A  spectrograph  was 
used  to  view  the  oxy  diborane  system  thtough  the  end  plate  of  a  single  pulse  shock  tube. 
Hydroxyl  and  several  boron  intermediates  were  identified,  qualitatively  verifying  both 
mechanism  and  equilibrium  calculations.  Ignition  induction  measurements  were  performed  up 
to  950° K  behind  the  reflected  shock  wave.  These  measurements  extend  previously  known 
data  by  two  orders  of  magnitude  in  the  induction  time.  Part  Two  of  this  report  describes  the 
CAL  Nonequilibrium  Normal  Shock  Wave  Program  which  was  extended  for  use  with  highly 
exothermic  branching  chain  reactions,  checked  out  with  known  hydrogen  oxyger  kinetics, 
and  reformulated  to  include  transport  and  radiative  energy  transfer  mechanisms. 
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INTRODUCTION  AND. SUMMARY 

Combustion  studies  reported  herein  follow  t lie  shock  tube  approach  toward  determination 
of  oxidation  mechanisms  and  reaction  kinetics,  applying  established  techniques  to  the  study  of 
pyrophoric  fuels.  These  studies  are  set  within  the  framework  of  a  long  range  AFOSR  project 
directed  toward  understanding  and  use  of  the  mutual  effects  of  chemical  reactions  on  collisional 
transport  properties  and  radiative  energy  transfer  in  flow  systems  of  high  power  density  heat 
release.  Results  of  the  Martin  Marietta  studies  include: 

1  A  branching  chain  mechanism  was  constructed  for  stoichiometric  oxygen/diborane 
(C^/BqH^)  systems  under  high  dilution  in  argon,  consistent  with  thermochemical, 
kinetic,  and  structural  data.  Included  in  fourteen  (forward)  steps  of  the  postulated 
mechanism  are  thermal  decomposition  of  diborane  and  production  of  hydroxyl  (OH) 
as  a  chain  carrier,  followed  by  the  plausibly  observable  oxides  and  hydrides  of  boron. 

2  A  single-pulse  chemical  shock  tube  was  used  to  study  the  high  temperature 
homogeneous  reactions  of  the  stoichiometric  02^2^  sYstc,n-  Ignition  induction 
measurements  were  obtained  for  the  pyrophoric  02/B9H()  system  in  the  temperature 
range  from  575  to  945°X  and  for  the  O2/CO  system  at  2,000°K.  These  measurements 
extend  the  previously  known  data  by  two  orders  of  magnitudes,  into  the  microsecond 
regime.  A  spectroscopic  emission  study  of  O2/B2H6  served  to  identify  BQ7  and  BO 
lines  in  the  spectrum,  taken  in  the  detonative  mode  at  temperatures  exceeding  945°K. 

3  Equilibrium  normal  shock  wave  calculations  were  performed,  supplementing  these 
emission  studies,  in  an  effort  to  establish  a  temperature  region  for  measurement  of 
intermediate  species  rate  constants  via  molecular  or  atomic  absorption 
spectrophotometry. 

4  An  OH  ultraviolet  absorption  experiment  was  designed  following  the  techniques  of 
References  1  and  2. 

5  Working  in  cooperation  with  Cornel!  University’s  Chemistry  Department,  the  Aero 
Combustion  Laboratory  at  the  Martin  Marietta  Orlando  Division  has  extended  and 
verified  the  use  of  Cornell  Aeronautical  Laboratory’s  Nonequilibrium  Normal 
Shock  Wave  Program  (Reference  3)  in  the  study  of  exothermic  systems  which 
include  branching  chain  mechanisms.  In  particular,  the  previously  measured 
appearance  and  disappearance  of  intermediate  species  occurring  during  the 
oxidation  of  hydrogen  from  about  1370  to  1540°K,  and  at  1  Atm,  (References  3 
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and  4)  was  successfully  simulated,  using  the  CAL  program  coded  for  the  Martin 
Marietta  OX*  6400  electroitie  eoniputer.  Martin  Marietta  calculated  ignition 
induction  times  are  compared  to  the  measurements  of  Reference  3. 

6  The  CAl.  program*  is  being  used  optionally  and  as  subroutine  to  a  Martin  Marietta 
developed  reaction  profile  model  which  incorporates  mutual  effects  of  chemical 
reactions  on  collisional  transport  properties  and  radiative  energy  transfer.  The 
transport  properties  subroutine  portion  of  the  Martin  Marietta  program  was 
satisfactorily  checked  out  for  various  multicomponent  mixtures.  A  fully  developed 
chemical  source  matrix  is  required  as  input  to  the  CAL  subroutine,  for  each 
reaction  mechanism  to  be  studied.  A  partially  known  or  postulated  mechanism  may 
be  described  parametrically  with  the  input  source  deck. 

Such  computer  methods  are  useful  to  the  chemical  kincticist  in  performing  diagnostics  on 
a  system  for  which  some  of  the  pre-exponential  factors,  activation  energies,  and  temperature 
exponents,  associated  with  unknown  specific  reaction  rates,  are  currently  being  determined  by 
measurement.  The  computer  program  and  various  sets  of  spcctrometric  measurements  are  used 
in  an  iterative  determination  of  mechanism  and  reaction  kinetics  consistent  with  minimum 
energy  and  detailed  balancing  principles.  This  technique  of  complementary  measurement  and 
numerical  experimentation  may  prove  superior  t<>  earlier  methods  used  in  attempts  to  elucidate 
mechanism  which,  although  yielding  results  consistent  with  measured  data,  did  not  serve  to 
guarantee  uniqueness. 

The  work  reported  herein  shows  how  steps  may  be  taken  to  isolate  a  physically  meaningful 
region  for  the  investigation  of  exothermic  reactions  using  a  shock  tube.  Brieuy,  this  is 
accomplished  by  establishing  the  onset  of  detonations  for  the  pressures  and  compositions, 
initially  set,  that  determine  shock  strengths  and  intermediate  compositions  to  be  monitored 
spectroscopically.  Though  the  present  work  did  not  proceed  to  the  next  step,  viz.,  measurement 
of  the  production  rates  and  concentrations  of  intermediates,  it  is  noted  that  techniques  in 
absorption  spectrophotometry  are  established  and  under  continuous  development.  It  is 
suggested  that  one  such  technique,  outlined  in  Appendix  A,  may  be  appropriate  in  extending 
the  present  work  to  the  proper  field  of  chemical  kinetics. 

in  the  mathematical  portion  of  this  investigation  it  is  shown  that  the  choice  of 
molecular  model  is  important  to  high  temperature  flow  field  applications.  Inasmuch  as  the 
molecular  dynamicist,  whether  he  be  engineer  or  kincticist.  now  requires  extensive  use  of 
computer  programs  to  follow  reaction  paths,  it  is  suggested  that  continued  effort  in  the 
development  and  use  of  programs,  such  as  the  one  herein  outlined,  will  yield  meaningful 
results. 


Part  One  of  this  report  presents  a  discussion  of  shock  tube  design  considerations  and 
preparation  of  experimental  facilities  directed  to  the  task  of  performing  shock  tube  kinetics. 

‘Herein  referred  to  as  the  CAl.  subroutine. 
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Although  no  rate  measurements  Inmc  been  made,  the  following  preliminary  tasks  were 
accomplished: 

1  Ignition  induction  times  were  measured,  extending  previously  known  data  by  two 
orders  ot  magnitude. 

2  Detonation  limits  were  determined  expcrimentallly . 

3  An  oxidation  mechanism  was  postulated  for  the  considered  pyrophoric  system. 

4  Equilibrium  compositions  were  calculated. 

5  Emission  spectra  were  photographed. 

6  Intermediate  chemical  species  were  identified  trom  the  spectra. 

7  An  OH  ultraviolet  absorption  experiment  was  designed  to  measure  specific  reaction 
rates. 

Part  Two  treats  the  general  mathematical  reaction  profile  model  under  development,  giving 
results  of  calculations  obtained  with  the  Martin  Marietta  modified  CAL  normal  shock  wave  and 
transport  properties  subroutines. 
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Shock  Tube  Investigations 


of  the  Diborane/Oxygen  System 
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I.  THE  SHOCK  TUBE  CONSIDERED 
A.  DESIGN 

The  shock  tube  was  designed  based  on  ideal  flow  conditions  using  equacions  from 
Reference  5.  The  shock  tube  is  divided  into  two  sections,  a  driver  section  filled  with  gas  of  low 
molecular  weight,  M4,  at  high  pressure,  p^,  and  a  driven  section  filled  with  an  inen  diluent  at 
low  pressure,  p^.  The  division  takes  place  at  a  diaphragm  burst  station  where  a  solenoid 
actuated  pin  ruptures  a  cellophane  diaphragm.  The  high  pressure  driver  gas  interface  on  low 
pressure  driven  gas  then  acts  as  a  piston  creating  a  shock  wave  which  moves  uniformly  down  the 
low  pressure  section  to  a  reflection  plane.  The  interface  (contact  surface)  between  the  driver 
and  driven  gases  follows  the  shock  wave  and  is  intercepted  by  the  reflected  wave.  Chemical 
reactions  among  the  reagents  in  the  inert  gas  take  place  during  the  time  between  reflection  of 
the  wave  and  its  interception  by  the  contact  surface,  the  observation  time.  Thus,  during  the 
observation  time,  spectroscopic  study  of  the  test  sample  is  carried  out  in  the  region  behind  the 
reflected  shock  wave  at  the  far  end  of  the  tube.  Analysis  of  the  shock  system  is  given  in 
Appendix  B. 

For  purposes  of  illustration  a  general  shock  tube  design  based  on  ideal  incident  shock  wave 
Mach  number,  obtained  at  infinite  diaphragm  burst  pressure  ratio  p4j  from  Equation  (32)  of 
Appendix  B,  is  presented  in  Figure  1.  The  solid  curves  represent  selected  cases  of  the  quadratic 
expression  for  the  incident  shock  Mach  number  variation  with  the  diaphr.  gm  burst  sound  speed 
ratio  for  parametric  combinations  of  74  and  7j,  the  ratios  of  specific  heats  for  the  driving  and 
driven  gases  respectively.  The  circled  point  marked  He/Ar  represents  the  asymptotic  shock  Mach 
number  attainable  when  helium  gas  at  very  high  pressure  is  used  to  drive  argon  at  very  low 
pressure,  the  gases  being  initially  at  room  temperature.  By  heating  the  driver  gas,  that  is,  by 
raising  temperature  T4,  one  may  expect  to  increase  the  shock  Mach  number  along  the  same 
curve.  The  locus  of  points  represented  by  the  dashed  line  shows  that  a  decrease  in  the  molecular 
weight  of  driver  gas  or  an  increase  in  molecular  weight  of  driven  gas  increases  the  shock  Mach 
number.  The  latter  technique  is  used  in  tailoring  the  gases  when  changes  in  molecular  weight  are 
introduced  by  mixing  test  reaction  samples  with  the  inert  diluent  (driven)  gas. 

The  choice  of  reaction  environment  fur  present  experiments  was  obtained  from  the  helium 
and  argon  combination,  a  finite  pressure  ratio  p^/p] ,  and  the  temperature  range  shown  in  Table 
I;  the  shock  Mach  number  does  not  exceed  a  nominal  value  of  3.6,  which  is  less  than  one  third 
of  the  asymptotic  value  shown  on  Figure  1,  A  graphic  representation  of  most  of  the  tabulated 
data  is  shown  on  Figure  2.  Note  that  the  total  shock  tube  length  xj+  x->  ^21  feet  is  nominally 
the  maximum  length  of  the  Orlando  Division’s  chemical  shock  tube. 
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Figure  1.  Asymptotic  Incident  Shock 
Mach  Number,  p^j/p^  -*■  oo 
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TABLE  I 


Shock  Tube  Parameters  at  T  j  s  298  15*K  for  Helium/ Arj»«*i 
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B.  EXPERIMENTAL  PROCEDURE 

Gas  dynamical  and  chemical  calibration  of  shock  wave  strengths  in  argon  and  with  argon 
diluted  oxidation  of  carbon  monoxide  was  conducted,  using  helium  as  the  driver  gas. 
Approximately  ideal  shock  tube  performance  was  obtained  experimentally  by  causing 
cellophane  and  alun  inuin  diaphragms,  inserted  between  the  high  and  low  pressure  sections  of 
the  shock  tube,  to  burst  at  pressure  t.dterences  approaching  the  material  yield  point.  The  effect 
of  small  differences  in  laboratory  environmental  conditions  on  ideal  pertormam  e  was  found 
negligible  from  both  experimental  and  theoretical  considerations.  The  experimental  ignition 
induction  period  was  identified  lor  the  oxidation  of  carbon  monoxide  using  Kistler  (>03  quartz 
crystal  transducers  to  measure  the  pressure  ise  due  to  combustion.  An  optical  system  using  an. 
RCA  031  A  photomultiplier  tube  was  designed  and  incorporat  d  into  the  shock  tube  viewing 
parts  to  measure  emission  outputs  during  the  ignition  induction  period.  Calculations  were 
performed  for  mixing  and  handling  the  oxy-diborane  system  in  the  experiments. 

1.  Instrumentation 

The  one-inch  inside  diameter,  single  puls  •  chemical  shock  tube  reactor,  its  associated 
vacuum  equipment,  electronics,  optical,  gas  handling  and  chromatographic  systems  which  have 
been  developed  and  used  during  the  contract  are  shown  in  Figure  3. 


5 


Figure  3.  Martin  Marietta 
Aero-Combustion  Laboratory- 


Equipment  used  in  this  laboratory  include  the  following: 

1  P.  E.  154D  Vapor  Phase  Chromatograph, 

2  Gas  mixing  and  vacuum  equipment, 

3  Shock  tube  and  36  liter  dump  tank,  with, 

a  Photomultiplier  mounted  on  the  reaction  head, 

b_  Control  panel  and  recording  equipment, 

4  Gertner  monochromator, 

5  Optical  and  mass  spectrometers. 

2.  Laboratory  Techniques 

The  gas  handling  system  is  used  to  mix  diborane  or  carbon  monoxide  and  oxygen  in 
stoichiometric  proportions,  using  argon  as  the  diluent.  When  mixing  pyrophoric  fuels  with 
oxygen,  the  reactants  are  mixed  separately  with  the  diluent  to  keep  the  reactants  below  their 
e  tpiosion  limits  at  room  temperature,  Referenced. 

a.  Events  In  The  Shock  Tube 

The  freshly  mi:  *cl  gases  are  allowed  to  stand  from  one  to  twelve  hours  before 
using,  one  hour  beinu  the  minimum  time  for  dilute  species  to  diffuse  a  distance  of  50  cm 
assuming  a  diffusion  coefficient  of  0.3  tm^  sec'^- 

After  allowing  the  mixture  to  diffuse,  the  low  p-  .e  section  of  the  tube  ischaiged 
with  a  mixture  sample  at  initial  pressure  pi .  Since  p}  may  vary  from  5  torr  to  5  psia  depending 
on  the  diaphragm  rupture  .<nd  incident  shock  Mach  number  requirements,  the  driven  section  is 
evacuated  to  about  1  micron,  using  a  rough  pump  when  these  conditions  require  that  p ]  S  5 
psia.  At  the  high  Mach  numbers  (pi  S  5  tofr),  the  2  inch  diffusion  pump  is  used  to  ensure  that 
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residual  air  molecules  in  t lie  driven  section  are  less  than  1  percent  of  the  number  density  of 
r  reactant  species.  An  argon  flush  is  used  after  each  run.  Pressures  in  both  the  high  pressure  (P4) 

and  low  pressure  sections  of  the  tube  are  measured  with  a  Bourdon  gauge  ca':,’rated  against  a 
standard  to  wLhin  1  percent.  Figure  4  shows  an  assembly  drawing  of  the  shock  tube  including: 

1  Dump  tank, 

2  Continuously  variable  high  pressure  p^  driver  section, 

3  Low  pressure  p|  driven  section  (the  reactor), 

4  Hi  vac  inlet  and  bypass  valves, 

5  Diaphragm  station  and  rupture  pin, 

_6_  Reactor  head  with  ball  valve  shutoff, 

7  Pressure  transducer,  product  sampling  and  optical  view  ports, 

Alternate  multiple  diaphragm  assembly. 


□  □ 


In  a  typical  run  sequence,  the  tube  is  evacuated  to  about  1  micron,  then  helium  is 
charged  into  the  high  pressure  section  to  about  1  kilotorr.  Samples  of  the  fuel  and  oxidizer 
selected  for  study  are  separately  premixed  with  an  inert  diluent,  argon.  In  the  non-pyrophoric 
experiments  with  CO,  the  reactants  arc  introduced  into  the  same  mixing  flask  with  the  diluent. 
The  mixtures  are  then  charged  into  the  tube’s  reaction  section  at  low  values  of  the  total 
pressure,  the  reactant’s  partial  pressures  being  such  as  to  minimize  their  collision  frequency.  The 
reaction  section  is  separated  from  the  high  pressure  driver  by  cellophane  or  aluminum 
diaphragms.  A  pin  is  caused  to  rupture  the  diaphragm  allowing  helium  at  high  pressure  to  enter 
the  reaction  section  creating  a  shock  wave  ahead  of  the  gas  interface.  The  shock  wave  excites 
the  gas  mixture,  elevating  its  temperature,  increasing  the  reactants’  collision  probability.  The 
shock  wave  then  reflects  from  the  end  plate  and  proceeds  into  the  excited  region,  elevating  the 
gas  temperature  a  second  time.  Thus,  homogeneous  exothermic  reactions  are  observed  in  the 
region  of  high  collision  probability  between  the  reflected  shock  wave  and  the  end  plate. 

Figure  5  shows  a  wave  diagram,  repiesenting  the  sequence  of  events  along  the  tube 
following  rupture  of  the  diaphragm.  Assuming  ideal  rupture  (instantaneous  removal),  the 
diaphragm  is  replaced  by  a  contact  surface,  gas  interface,  which  accelerates  into  the  low  pressure 
section.  During  its  short  period  of  acceleration,  small  disturbances  project  forward  from  the 
contact  surface  (at  the  acoustical  speed  in  Region  1).  These  small  disturbances  coalesce  into  a 
shock  wave  after  which  shock  wave  and  contact  surface  both  move  forward  at  uniform  velocity. 


Prtuurt 
Distribution 
At  Ttmaa 
t»tj  and  t*ta 

Aftar 

Diaphragm 

Ruptura 


Figure  5.  Chemical  Shock  Tube  Wave  Diagram 
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Coincident  with  formation  of  the  shock  wave,  a  rarefaction  (expansion)  wave  moves 
from  the  decelerating  contact  surface  toward  the  driver  section  end  wall,  reflects  from  the  end 
wall,  and  is  accelerated  toward  the  reaction  section,  cooling  the  gas  in  Region  3.  Provided  the 
tube  operating  conditions  are  properly  selected,  the  head  of  the  reflected  expansion  wave 
terminates  on  the  contact  surface-reflected  shock  wave  intersection.  In  this  case  the  shock  wave 
is  partly  transmitted  into  t lie  cold  driver  gas.  Region  3,  and  partly  reflected  into  Region  5  as  an 
expansion  wave  which  cools  the  twice  processed  gas,  effectively  quenching  the  reaction.  The 
transmitted  shock  wave  is  attenuated  by  the  open  dump  tank,  thereby  preventing  further 
heading  of  the  gas  sample.  The  reaction  quency  rate  stated  in  the  literature  is  1  ()^° K/sec  over  an 
average  dwell  time  of  a  few  hundred  microseconds.  Because  of  this  effective  quench  (Reference 
7)  the  tube  is  called  a  single  pulse  chemical  shock  tube. 

b.  Non-Idealities 

This  technique  of  observing  exothermic  reactions  in  the  relatively  quiescent  region 
downstream  of  the  reflected  shock  wave  (Region  5)  is  being  used  extensively  (Reference  8 
through  14).  The  princip  d  advantage  considered  is  that  measurements  are  facilitated  at  high 
temperatures  not  available  with  the  incident  shock  wave  technique  (measurements  made  in 
Region  2)  because  detonations  likely  incurred  thereby  tend  to  complicate  the  measurement. 
However,  caution  should  be  exercized  when  making  rate  measurements  with  a  single  pulse  tube. 
This  is  because  of  certain  flow  non  idealities  as  follows  (Reference  15): 

J_  Diaphragm  rupture  causes  large  initial  flow  perturbations  that  could  be  amplified 
if  there  is  large  exothermicity  in  Region  2. 

2_  Boundary  layer  build-up  at  the  reflection  end  plate  tends  to  cool  large  portions 

of  the  test  samples  when  these  are  observed  far  from  the  end  plate. 

3  Slow  cooling  by  means  of  the  centered  rarefaction  wave  may  cause  systematic 

errors  in  the  data  because  the  effective  reaction  test  time  must  be  integrated  over 
the  cooling  time  (Reference  16  and  17). 

Caution  exercized  in  the  present  work  circumvents  many  of  these  nonideal  flow 
problems.  The  first  problem  is  alleviated  in  proportion  to  the  amount  of  diluent  used:  96 
percent  Argon  for  the  stoichiometric  oxy-diboranc  mixture  studied.  Evidence  that  detonations 
did  not  occur  in  Region  2  while  measuring  ignition  induction  times  is  given  below.  Further, 
according  to  the  criterion  given  in  References  15  and  18,  exothermicity  below  that  required  for 
detonations  in  Region  2,  but  which  could  otherwise  cause  flow  perturbations  sufficient  to 
accelerate  the  incident  shock  wave,  is  generally  avoided  by  ensuring  that  the  incident  shock 
wave  velocity  uj  is  greater  than  the  Chapman-Jouguct  (CJ)  detonation  velocity  uqj  of  the  test 
mixture.  The  method  outlined  in  Appendix  B  was  used  to  compute  the  equilibrium 
concentrations  of  species  in  Regions  2  and  5  under  test  conditions.  The  species  mole  fractions 
other  than  the  initial  reagents  and  diluent  were  found  to  be  negligable  in  Region  2  under  those 
conditions.  The  flow  field  is  therefore  frozen  across  the  incident  shock  wave,  with  7]  =  1 .634.  It 
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follows  that  u  |  ’U(/|  when  the  incident  slun  k  Mach  number  is  greater  than  the  square  root  of 
the  temperature  ratio  across  the  slun  k,  i.  e..  when  M|  >  y/  I 1  |  ■  I  his  was  indeed  the  ease  for 
diborane  ignition  induction  measurements  presented  below. 

The  second  problem  is  alleviated  according  to  Mark’s  criterion.  Reference  1  9,  viz.  that 
the  effective  specific  iieat  ratio  be  7|  >  1.5  so  that  the  extent  of  test  sample  cooling  due  to 
boundary  later  interaction  with  the  bifurcated  shock  wave  is  diminished.  This  is  also 
accomplished  by  placing  the  observation  port  close  to  the  end  wall  of  the  tube.  In  the  tube 
used,  the  ratio  of  port  to-end  wall  distance  over  the  tube  diameter  is  0.74. 

The  third  problem  is  alleviated  by  using  the  single  pulse  shock  tube  modification 
introduced  by  Klepeis,  Reference  20,  viz.  the  use  of  a  dump  tank  designed  to  increase  the 
cooling  rate  by  factors  of  from  5  to  10,  Reference  17. 

3.  Performance 

a.  Transducers  and  Shock  Velocity 

The  shock  wave  velocity  has  been  measured  variously  with  Martin  Marietta  desig  ed 
barium  titanate  and  quartz  piezoelectric  pressure  sensors  and  by  Kistler  701.  606L,  and  603 
pairs  of  quartz  pressure  transducers.  The  acceleration  compensated  Kistler  603’s  give  the 
sharpest  response,  largest  signal-to-noise  ratio  and  most  ideal  (response  to  a  step  input)  damping 
characteristics  of  all  transducers  tested;  see  the  lower  beam  traces  of  Figure  6  c  and  d. 
Transducer  spacing  was  designed  to  obtain  average  velocities  at  various  tube  stations  located 
between  the  reflection  end  plate  and  the  diaphragm.  A  delay  signal  is  triggered  by  a  barium 
titanate  transducer  mounted  closest  to  the  origin.  This  sensor  was  chosen  for  the  trigger  because 
of  its  superior  sensitivity.  The  delay  signal  is  swept  at  a  speed  chosen  so  as  to  display  a 
maximum  amount  of  information  on  the  recording  trace.  Passage  of  the  shock  wave  is  observed 
on  the  traces  of  Figure  6.  The  two  upper  oscilloscope  records  (a  and  b)  correspond  to  an 
initially  set  pressure  ratio  pq/p]  =  40  (idea!  incident  shock  Mach  number,  M  =  3.075).  Record  a 
is  a  helium/argon  run  with  Kistler  606L’s  sensing  the  shock  wave  between  transducer  stations 
separated  10  inches.  Note  the  short  time  lapse  between  the  second  and  third  pressure  rise. 
Approximately  1/4  of  the  lapse  corresponds  to  the  time  between  arrival  of  the  shock  wave  at 
the  second  transducer  (the  second  rise)  and  the  time  of  reflection  of  the  incident  shock  wave 
from  the  end  plate  of  the  tube.  Since  the  reflected  shock  wave  moves  at  about  one-half  the 
speed  of  the  incident  shock,  the  third  time  lapse  (between  the  3rd  and  4th  rise)  lasts  about 
twice  as  long  as  the  first.  This  is  better  observed  on  the  three  lower  records  which  correspond  to 
an  initially  set  pressure  ratio  P4/ p ^  =  20  (M  =  2.55). 

b.  Ignition  Induction  Check  Runs 

Notice  the  flat  (but  noisy)  structure  of  the  third  plateau  on  record  6a.  Record  6b 
corresponds  to  a  30  percent  stoichiometric  mixture  of  ()7  and  CO  diluted  in  argon. 
Similar  tests  have  been  reported  by  K.G.P.  Sultzmann  (Reference  21).  Evidently,  the  5  percent 
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(a)  Upper  Trace:  1000  u  sec/cm  (b)  Upper  Trace:  1000  u  sec/cm 

Lower  'Trace:  100  u  sec/cm  lower  Trace:  100  u  sec/cm 

Iielium/Argon  tare  run  Helium  Argon  {}0%  CO-O.  ) 

Pp  =  TO  mm  Hg  Px  =  20  mm  Hg  2 


(c  )  Upper  Trace  -  i  hotomu.lt  iplier  output 
Lower  Trace  -  Transducer  Out;ut 
sweep  rate  -  /0  u  sec/ cm 
Helium/ argon  taro  run 
la  a  lOg.L  mm  Hg 


— 


(e) 


Lower  Trace  -  Transducer  Cutout 
Lweep  rate  -  50  u  sec/cm 
Pp  =  1Q3.4  mm  Hg 

Helium/„rgon  (3Ct©  CO  -  0^) 


Upper  Trace  -  Photomultiplier  Output 
Lower  Trace  -  Transducer  Output 
sweep  rate  -  200  u  sec/cm 
h'eliuVArgon  (go,©  CO  -  0  ) 

Pi  =  10 4  mm  Hg 


Figure  (i.  Carbon  Monoxide  Flame  Experiments 
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change  in  specific  heats  lias  hail  very  little  or  no  effect  on  the  incident  .-.hock  velocity  or  on  the 
level  of  the  third  plateau,  compared  to  the  tare  run.  However,  the  effect  of  combustion  is 
evident  on  the  third  plateau  where  a  superposed  pressure  due  to  the  oxidation  of  CO  can  be 
observed  to  span  about  250  microseconds.  The  ideal  temperature,  T^,  in  the  reaction  zone, 
behind  the  reflected  shock  wave,  is  2270°K  tor  P4/ p -j  -  40,  which  e'cceds  the  auto-ignition 
temperature  of  CO  (924° K)  by  a  factor  of  about  2'/>,  The  actual  temperature  is  closer  to 
2000° K,  determined  from  the  Martin  Marietta  calibration  curve  (discussed  below)  and  from  the 
theoretical  data  on  Figure  2.  Hence,  ignition  induction  commences  at  the  shock  front  and  lags 
about  250  microseconds  piior  to  maximum  burning.  This  is  in  agreement  with  the  results  oi 
Sultzmann.  Note  that  p^  =  20  torr  for  these  runs  (rather  than  103.4  torr),  accounting  for  a 
reduced  collision  probability,  hence  no  observable  pressure  rise  due  to  auto  ignition  temperature 
behind  the  incident  shock. 


c.  Optics 

An  attempt  to  verify  the  ignition  induction  period,  by  observing  the  total  CO  flame 
emission,  was  implemented  by  designing  an  optical  system  (Figure  7)  with  an  RCA 
931 A  photomultiplier  tube  located  90  degrees  azimuthaily  from  the  end  transducer  in  order  to 
view  light  emitted  from  the  slug  of  gas  immediately  downstream  of  the  .eflected  shock  wave. 
The  sequence  of  oscilloscope  records,  (Figure  (>,  c,  d,  and  c)  shows  a  hclium/argon  tare  run,  a 
helium/ 30  percent  CO/0->  stoichiometric  mixture  diluted  in  argon,  both  swept  at  50  jUscc/cm, 
and  a  repeat  of  the  oxidation  run  swept  at  200  pscc/cm  in  order  to  display  the  entire  history  of 
events,  all  for  pq/pj  =  20.  As  in  the  case  of  the  pq/p]  -  40  runs,  burning  raises  the  level  of  the 
third  and  fourth  plateaus  whereas  no  change  is  observed  nor  predictable  on  the  level  of  the  first 
and  second  plateaus,  the  temperature  in  this  case  being  about  130°K  less  than  the  auto-ignition 
temperature.  However,  in  this  case  the  initially  set  total  charge  pressure  was  about  5  times  that 
of  the  earlier  pq/p^  -  40  runs,  (Figure  6a,  h).  That  is  p^  =  103.4  torr  so  that  the  combined 
temperatute  and  pressure  effect  can  be  observed  in  the  initial  photoemission  rise  (negative 
polarity,  upper  traces)  over  a  period  of  about  10  microseconds.  Evidently,  the  ignition 
induction  period  is  complete  by  the  time  that  the  shock  wave  lias  arrived  at  the  end  transducer 
after  reflection  from  the  end  plate. 


d.  Shock  Tube  Calibration 

Approximately  ideal  air/air  shock  tube  performance,  shown  in  Figure  8  was 
accomplished  experimentally  after  obtaining  near  ideal  diaphragm  rupture  (pctaling)  using 
cellophane  and  scored  aluminum  foil  diaphragms.  This  technique  requires  that  pq/pj  be 
maintained  at  diaphragm  pressure  differences  approaching  the  material  yield  point.  Thus,  the 
diaphragm  ruptures  with  a  minimum  of  material  particles  being  lost  to  the  stream  and  with  a 
corresponding  minimum  energy  transfer  from  the  stream  to  the  diaphragm  particles.  The 
transducer  separation  for  those  data  was  17.80  cm. 


<L 


*  Grind  Spherical  Radius  -  0.56  in 
Aluminum  Mirror  Surface  Polished 


Figure  7.  Optical  Kmission  Sensor 


Attempts  have  been  made  to  eliminate  the  departure  from  ideal  for  the  Hc/Ar 
calibration  by  running  experiments  at  various  transducer  separations  to  increase  the  resolution 
(obtain  an  optimum  ratio  of  transducer  sensing  element  diameter  to  transducer  spacing).  The 
best  resolution  was  obtained  with  the  smallest  transducer  spacing  tried.  Consequently,  the  tare 
runs  for  a  transducer  spacing  of  6.68  cm  was  used  to  define  the  standard  calibration  curve.  This 
curve  is  shown  in  Figure  9  where  the  present  experiments  are  compared  to  both  theory  and 
results  of  previous  investigations. 


c.  Sources  of  Error,  Data  Scatter,  and  Uncertainties 

The  principal  sources  of  error  contributing  to  the  departure  of  experimental  shock 
tube  performance  from  the  ideal  are: 

1  Sensor  resolution 

2  Non-ideal  diaphragm  rupture 

3  Boundary-layer  effects 

4  Laboratory  environmental  variations 


Attempts  to  eliminate  the  first  two  sources  of  error  listed  have  been  discussed  above. 
A  correction,  in  the  case  of  non-ideal  rupture,  is  suggested  but  has  not  been  attempted  at  this 
writing.  It  is  suggested  that  the  diaphragm  may  be  weighed  on  an  analytical  balance  both  before 
and  after  a  run.  The  energy  balance  required  to  accelerate  diaphragm  residue  down  the  tube 
may  then  be  used  to  correct  for  shock  wave  attenuation. 

Shock  wave  attenuation  due  to  the  boundary  layer  is  another  factor  which  can  be 
accounted  for  in  order  to  correct  the  departure. 

The  data  shown  on  Figure  8  indicates  good  repeatability,  considering  that  each  of  the 
single  data  points  plotted  at  a  given  pq/p^  represents  two  indistinguishable  runs. 
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Figure  9.  Helium/Argon  Performance  Calibration 
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II.  TECHNICAL  APPROACH 
A.  SELECTION  OF  FUELS 

The  selection  of  high  energy  fuels  used  in  the  experimental  studies  was  predicated  on 
existing  thcrmochemical  and  kinetics  data,  commercial  availability,  and  ease  of  synthesis. 
Prior  literature  survey,  dv.  cribed  in  Reference  22,  eliminated  pentaborane  and  aluminum 
borohydride  from  1  nmediate  practical  consideration  because  of  their  high  costs.  A  considerable 
amount  of  research  has  been  conducted  on  the  boron  hydrides,  their  stable  adducts,  and 
substitution  alkyl  compoui.ds,  (Refeiomes  23  through  27).  In  Reference  23,  for  example,  the 
aim  of  research  has  been  to  obtain  a  fundamental  understanding  of  combustion  and 
extinguishment  for  the  d.borrne-oxygen  (B2H^/0'>)  and  for  the  pentaoorane-oxygen 
reactions.  Ignition  delay  times  for  82H^  anc^  BcjHq/oxygen  mixtures  diluted  in 
argon  were  measured  using  a  single  pulse  shock  tune  adapted  to  the  tailored  interface  technique 
reported  in  Reference  28.  Shock  lube  ignition  of  the  B^Hq/C^/Ar  mixtures  required  slightly 
higher  temperatures  than  those  coi  raining  diboranc,  although  static  bulb  experiments 
1  Reference  23)  indicate  that  pentaborane  is  more  flammable  than  diborane. 

i  cr  both  compounds,  test  results  show  that  ignition  delay  times  are  independent  of  the 
oxygen  concentration  throughout  an  equivalence  ratk  range  from  0.5  to  2.  This  information  is 
based  on  an  equivalence  ntio  of  1  corresponding  to  a  1:6  pentaborane-oxygen  ratio  which 
could  react  completely  according  to 

2  B5H9  +  12  02— ►  5B203  +  9  H20: 
and  corresponding  to  a  1 :3  diboranc  ratio  for  the  overall  reaction 

B2H6  '  3  02—  B2(>3  +  3  H20  . 

Because  of  the  oxygen  concentration  independence  of  the  ignition  induction  times  and  since  the 
activation  energies  are  close  to  those  for  pyrolysis,  it  was  concluded  in  Reference  23  that 
pyrolysis  is  the  rate  controlling  step  in  ignitic" 

Figure  10  presents  thcrmochemical  data  in  favor  of  the  selection  of  B2H^  as  a  candidate 
for  oxidation  kinetics  studies.  Figure  10  shows  the  heat  liberated  in  combustion  of  1 :3  mixtures 
of  B2H^/02  over  the  temperature  range  300-2300°K,  compared  to  the  heat  liberated 
in  combustion  of  a  2:1  mixture  cf  H->/t  < .  (data  taken  from  the  JANAF  tables,  Reference  29). 
Thus,  in  applications  where  heat  liberation  per  pound  of  total  reactants  is  important,  the 
reaction  yielding  boron  oxide  in  the  crystalline  state  produces  1300  Btu  per  pound  more  than 
the  H2/02  reaction  at  300° K  whereas  the  excess  at  2200° K  drops  to  570  Btu  per  pound  for 
B203  in  the  liquid  state.  Applications  to  supersonic  combustion  seem  desirable. 
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of  Reactants 


Clearly,  the  best  way  to  proceed  with  a  diagnostic  method  of  investigation  which  one 
hopes  to  apply  to  the  understanding  of  mechanisms  and  kinetics  of  a  chemical  system  is  to  sec 
how  that  method  has  been  applied  to  the  understanding  of  other  systems.  Literature  search 
and  problem  solving  began  with  the  relatively  well  known  hydrogen/oxygen  (H2/02)  system, 
including  brief  digressions  with  the  better  known  hydrogen/halide  (H2/X2)  and  the  nitrous 
oxide/hydrogen  (N7O/H9)  systems,  all  of  which  are  branching  chain  mechanisms  at 
temperatures  above  600°K.  The  Martin  Marietta  reaction  profile  model  computer  program  was 
used  to  try  an  H7/O7  reaction  test  case.  Reaction  profile  models  are  now  part  of  the  molecular 
dynamicist’s  diagnostic  tool  kit,  both  in  industry  and  at  universities.  The  Orlando  Division  of 
Martin  Marietta  worked  with  Cornell  University  during  the  simultaneous  development  of  very 
similar  reaction  profile  models,  thereby  gaining  insight  to  the  problems  that  will  arise  when  the 
Orlando  Division  proceeds  to  more  complicated  systems  such  as  H2/air.  B2Hg/02  and  B2H^/air 
for  computer  diagnostics  or  engineering  applications.  The  details  of  this  work  will  he  given  in 
P;irt  Two  of  this  report. 


Experimentally,  the  B^H^/O^  system  is  more  complex  than  those  aformentioned  because 
element  number  five,  boron,  has  three  electrons  to  share  in  combination  with  other 
elements,  in  contrast  to  more  stable  elements  such  as  its  neighbor,  carbon.  Diborane  has  six 
pairs  of  electrons  to  distribute  among  nine  bond  positions,  that  is,  2/3  of  a  single  (pair)  bond 
per  bond.  The  molecule  is  therefore  said  to  be  electron  deficient  (on  the  average,  1/3  electron 
pair).  This  is  'he  fundamental  reason  for  the  reactivity  and,  consequently,  sometimes  desirable 
quality  of  diborane  as  a  high  energy  pyrophoric  fuel.  Also,  this  system  has  ten  individual  atoms 
(B-jH^/O^  to  permute  within  the  limits  of  certain  allowed  combinations.  These  limits  may 
change  in  a  way  that  is  yet  to  be  discovered.  In  contrast,  after  a  decade  of  modern  techniques 
including  the  use  of  shock  tubes  and  reactor  profile  analysis,  the  H2/02  system  is  now 
essentially  fixed  at  eleven  steps  in  its  high  temperature  and  homogeneous  branching  chain 
mechanism,  with  only  four  individual  atoms  to  permute! 


B.  HYPOTHESIZED  B2H6/02  MECHANISMS 

It  is  generally  conceded  in  the  liteiature  that  B2H^/()2  will  exhibit  a  branched  chain 
mechanism,  but  even  the  early  hypotheses  (e.g.  Roth's  thesis.  Reference  30)  have  not  been 
clearly  verified  beyond  the  first  step,  viz.,  1/2  B2H^  +  02-«-L BHj  +  02~- *■ 

Roth  and  Bauer  (References  31  and  32)  postulated  a  mechanism  for  the  B2H6/02  reaction 
(Reference  15)  based  on  partial  analogy  with  the  H?/02  at  the  second  explosion  limit,  namely 

B2H6^-2  BH3 
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bh3  +  o2— bh2oh  +  o 

O  +  B2H6-"BH2OH  +  bh3 
BH^OH  +  B2H6—  BH3  +  B2H5OH 
BH3  +  °2  +  M-HB02  +  H2  +  M 
BH3  +  B2H^-*-B3H2  +  H-) 

B3H7-  polymerization 
B2H5  OH— ' ►further  oxidation 

where  M  in  the  fifth  reaction  step  is  any  third  body.  The  mechanism  was  further  modified  by 
Goldstein  (Reference  33)  and  has  been  studied  by  the  group  at  Rensselaer  under  Bauer  and 
Wiberly  (Reference  34).  Compared  to  the  oxidation  of  hydrogen,  the  oxidation  of  diborane 
requires  considerable  study  before  its  mechanism  is  understood  as  fully  as  H2/02. 


Recently,  the  shock  tube  work  of  Skinner  ct  al  (Reference  23)  opened  a  new  line  of 
research  for  understanding  the  mechanism  of  flame  inhibition  in  both  H2/02  and  B2H^/02 
reactions.  The  investigations  of  Koski  at  Johns  Hopkins  and  Porter  at  Cornell  have  yielded  low 
temperature  mass  spectra  of  BH3  and  other  intermediate  species  in  the  B2H^/02  system 
(References  35  and  36). 


Gobbett  and  Linnett  (Reference  37)  have  studied  the  B2H6/02  system  in  a  conventional 
reaction  vessel  at  temperatures  between  393  and  483° K,  while  Jolly  and  Schmitt  (Reference  38) 
have  studied  the  reacion  of  diborane  with  hydroxide  in  aqueous  solutions,  lending  credence  to 
the  existence  of  BH(OH)2  as  an  intermediate. 


Based  on  stoichiometric  data  derived  in  these  experiments,  and  from  the  available 
literature  on  the  B^H^/O^  system,  the  following  reaction  mechanism  was  postulated, 
consistent  with  thcrmochemical  data  shown  on  Table  11: 
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(AH)  298°K’  kcal/mole 


B9H 


2n6^>  2BH, 


h2/o2 


t)2  +  BH3— >  H3BOO«=±H2B(X)H 

k 


Ar  +  BH9O 


2  BO  +  H9  +  Ar 
— >  e- 


Ar  +  OH  +  BO 


OH  +  BHO9 


\3_^  OBOH  +  Ar 


k , 


BO  2  +  H20 


(4  Center)  02  +  BH9 


OH  +  BH 3  k5  BH2  +  H20 
k, 


‘6  ^  BHO  +  OH 


02  +  BHO  k7  ^  BH02  +  O 
k 


h2  +  B02 
2  BH02 
O  +  H2 

Ar  +  2H 


8  ^  BH02  +  H 


c9  B9O3  +  H20 


‘10  OH  +  H 


11. 


H0  +  Ar 


Net:  3()2  +  B2H6 — >  B2Oj  +  3H20  (-386) 

Observables  /OH;  BO;  BC)o:  BH-^ 


:1 - »  H2BO  +  OH 

-154 

-10 

-45 

-38 

-75 

-7 


-104 


-442.0 


-387.0 


+36.0 


+  6.0 


+  1 1.0 

+  2 


+55.0 


Note: 


Prof.  S.  H.  Bauer  suggests  that  B02  +  B02H  is  more  likely  for  the  production  of 

b2o3. 


TAB1.F  II 


Thcrmochemical  Data.  H°j-‘ 


"K 

Molecule 

298 

1000 

1500 

B 

133 

133 

132 

boron 

BH 

106 

105 

104 

boron  monohydride 

BHO 

-20+20 

-21 

-21.5 

boron  oxide  hydride 

BHO-, 

134 

- 1  36 

-137 

metaboric  acid 

BFb 

48±  1  5 

47 

46 

boron  dihydride 

BH 

B?H,0 

1  1 4±  1 . 

115 

115 

boron  dihydroxide 

| 

bh3 

25. 5±  1 0 

22.6 

21.5 

boron  trihydride  | 

BH  3O 

-70.0+1 

(J  Phys.  Chem.  68, 

monohydroxyborane  1 

bh3o2 

-153.8 

27,  32,  64) 

dihydroxyborane 

BH  30  3 

-237 

240 

-241 

boric  acid 

B2 

1  95±6 

194 

192 

boron 

B-jH'jO  4 

-200.4±3.5 

(j.  Phys.  Chem.  68,  3' 

64,  '65) 

B2H6 

b2h3o 

9.8±4 

5.5 

5.7 

diborane 

B-?0? 

- 1  ()9±2 

-109.5 

-1 10.5 

dimeric  boron  monoxide 

b2o3 

-199.14 

200 

-201 

boron  oxide 

B  4H  3t  >  4 

-291 ±10 

293.6 

-293.5 

boroxin 

b3h3o6 

-543±3 

-542 

-540 

trimeric  metaboric  acid 

BO 

1 1  ±  1 0 

9.88 

8.76 

boron  monoxide  ! 

bo2 

-75.27 

-75.3 

-75.6 

boron  dioxide 

b2h4o4 

-307 

-309 

-308 

.  — 

boron  dihydroxide  (dimeric) 

*in  kcal/inole,  JANAF 


Reaction  steps  in  the  chain  include  the  thermal  decomposition  ol  B^H^,  witli  dissociation 
energy  of  1.62  cV  (Reference  39)  to  give  two  moles  of  BH3.  Oxygen  combines  with  13H3 
producing  one  of  several  isomers,  H3BOO^X  H->BOOH  at  high  temperatures, or  H B(OH)-j  at  low 
temperatures.  Following  Berl  (Reference  -40)  this  step  may  contribute  to  the  production  of 
hydroxyl  radical  and  H->BO.  Thermochemical  data  for  the  latter  molecule  is  not  yet  available.  It 
was,  therefore,  necessary  to  add  the  reaction  steps  leading  to  the  production  of  boron 
monoxide,  deleting  the  isomerization,  to  obtain  the  heat  of  reaction  of  6.0  kcal/mole.  The  total 
heat  liberated  in  the  sequence  of  reaction  steps  summed  to  within  0.26  percent  of  the  heat 
liberated  in  the  net  reaction. 


The  only  li  i*^l  i  t  cm  peril  tu  iv  rate  obtained  at  th  is  writing  is  the  one  for  pyrolysis.  Reference 
23.  Subsequent  reaction  steps  which  contribute  to  ignition  induction  may  he  determined 
tin  •uugn  a  search  tor  readily  observable  free  radicals  such  as  OH  or  boron  monoxide.  Thus,  a 
speetrospopic  investigation  of  the  ignition  induction  mechanism  was  initiated. 
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III.  SHOCK  INDUCED  IGNITION  OF  DIBORANE 


The  shock  tube  used  to  measure  ignition  induction  time  is  described  in  Section  I.  A 
schematic  of  the  shock  tube,  instrumentation,  and  recording  setup  is  shown  on  Figure  1 1. 


K  is  tier 
Charge  Ampl 


Figure  12  illustrates  the  region  of  observation  on  a  time-distance  wave  diagram, 
representative  of  events  occurring  axially  in  the  shock  tube  following  the  time  of 
diaphragm  burst. 

Measuring  station  (1)  is  defined  at  1.88  cm  from  the  shock  reflection  end  plate  of  the 
driven  tube.  Overall  length  of  the  driven  tube  is  40  inches,  from  the  diaphragm  station  to 
the  end  plate,  for  all  subsequent  experimental  data  presented.  The  driver  tube  length  is 
continuously  varied  (Reference  7)  to  adjust  the  wave  intersect  position  for  optimum  gas 
dynamic  quench  characteristics  (i.c.,  cooling  rates  of  lO^K/scc).  Thus,  the  total  shock  tube 
length  can  be  varied  from  about  4  feet  2  inches  at  the  high-diaphragm  burst-pressure  ratio,  P41 
-  100,  to  about  7  1/2  feet  at  the  low  P4  |  =5.  Helium  was  used  as  the  driver  gas,  to  produce 
shock  waves  in  argon. 
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Figure  12.  Wave  Diagram  -  Flame  Schematic 


An  RCA  931 A  photomultiplier  (PM'  is  mounted  at  an  observation  port  centered  on 
measuring  station  (1).  A  Kistler  603  pressure  transducer  *is  also  mounted  at  station  (1), 
perpendicular  to  the  observation  port,  so  that  both  light  and  pressure  intensities  arc  observed 
simultaneously  from  the  same  homogeneous  gas  sample  be!  d  the  reflected  shock  wave 

The  flame-like  schematic  represents  spatial  extent  of  the  intensity  observed  as  a  function 
of  time,  superimposed  on  the  wave  diagram.  (Such  flames  have  been  observed  in  shock 
tubes  using  Schlieren  photography.) 
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A.  IGNITION  INDUCTION  MFASlIk F.  ME  NTS 


Definition  of  ignition  induction  time,  t|,  for  the  present  measurements  is  shown  as  the 
incremental  time  from  arrival  of  the  reflected  shock  wave  at  the  measuring  port  to  the  time  of 
maximum  light  intensity  at  the  measuring  port  (Figure  12). 


Oscillogram  traces  from  the  experiments  are  shown  in  Figure  13  for  the  sequence  of 
diaphragm  burst  pressure  ratios  P.j  ^  =  5.  5.5,  6.  0.5,  7.  8,  8.9,  and  11.42,  representing  ideal 
reactor  bath  temperatures  in  the  range  750  to  1  175~K  behind  the  reflected  shock  wave.  The 
actual  temperature  range,  after  bulk  correction  fur  losses  due  to  boundary  layer  and  inefficient 
diaphragm  burst,  is  believed  to  be  575  to  945”  K  ±50  K.  The  tolerance  is  obtained  from 
literature  values  for  che  attenuation  of  shock  waves  in  shock  tubes  of  L/D  approximately  that 
used  in  this  experiment.  Data  reduction  for  this  series  of  runs  is  displayed  on  Table  III. 


Figure  13.  Oscillogram  Traces.  Ignition  Induction  Measurements  of 
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The  incident  shock  wave  velocity  is  measured  as  the  difference  between  the  first  two 
pressure  steps  on  the  lower  trace.  The  end  of  the  wiggle  between  the  second  and  hire  ste  "'s 
indicates  the  acoustical  effer*  cf  .hock  wave  reflection  at  the  end  plate.  The  third  and  fourth 
steps  measure  the  reflected  shock  wave  velocity.  Pressure  oscillations  after  this  time  are  due  to 
combustion.  Note  correlations  with  the  (upper)  photomultiplier  trace. 

To  obtain  measurements  under  approximately  ideal  conditions,  the  stoichiometric 
ru  rtants  =  1:3)  are  premixed  at  high  dilution  ratios  in  argon.  Initial  conditions,  i.c., 

the  initially  set  pressures  in  the  driver  and  driven  (realtor)  sections  of  tne  shock  tube,  arc  then 
set  according  to  the  theory  governing  propagation  of  ideal  shock  waves  in  the  inert  diluent 
(Appendix  B).  Both  measurement  and  theory  show  there  is  negligible  error  (Figures  14  and  15) 
in  this  procedure  provided  that  diluent  concentrations  are  sufficiently  la*ge  and  that  the 
incident  shock  wave  is  sufficiently  weak  to  preclude  the  onset  of  detonation  modes.  It  follows 
that  the  equilibrium  temperature  To  downstream  of  the  shock  wave  (Region  2  in  shock  wave 
nomenclature)  is  too  small  to  provide  significant  interparticle  collision  frequency  for  ignition  of 
the  highly  diluted  fuel/oxygen  gas  mixture.  This  is  verified  experimentally  with  the  PM  which 
views  the  incident  shock  wave  in  the  reaction  section  of  the  tube.  Since  no  light  is  detected  in 
the  region  of  sensitivity  of  the  PM,  centered  about  0.42  micron,  until  after  shock  reflection 
from  the  end  plate  and  since  there  is  no  measurable  difference  between  the  incident  sho  k  wave 
velocity  in  the  mixture  or  in  the  pure  diluent  for  the  same  initial  conditions,  it  follows  that 
ignition  occurs  only  after  shock  reflection  for  those  initial  conditions.  (The  observed  effects 
here  are  entirely  analogous  to  the  observations  with  CO  ignition.) 


Incident  Shock  Velocity  ^nim/psec 


Figure  14.  Shock  Waves  in  Pure  Argon,  Initial  Temperature  =  2P8°K 


Figure  13.  Shock  Waves  in  Ar  (96  Percent)  +  O2  (3  Percent)  +B^Hg 
(1  Percent);  Initial  Temperature  =  298°K;  Pressure  =  500  mmHg 


Under  these  conditions  it  may  he  assumed  that  the  initial  mixture  mole  fractions  are 
unchanged  in  traversing  the  shock  wave  into  Region  2.  Figure  1 5  shows  the  equilibrium 
chemistry  temperature  calculations  (Appendix  B)  which  assume  first,  that  flow  is  chemically 
frozen  across  the  incident  shock  wave  and,  second,  that  equilibrium  concentrations  of 
intermediate  species  postulated  in  the  chain  mechanism  on  page  occur  downstream  of  the 
reflected  shock  wave.  (Note  that  initial  pressure  variation  does  not  significantly  change  T5.) 
Martin  Marietta  experimental  data  is  plotted  on  Figure  1  5.  Comparison  of  these  data  show  that 
the  reported  experimental  temperatures  lie  above  the  equilibrium  values  by  about  40° K  at 
the  highest  measured  velocity.  This  is  not  surprising  in  view  of  the  fact  that  the  experimental 
temperatures  are  deduced  from  ideal  shock  calculations  for  the  pure  diluent,  Figure  14. 
Converting  the  measured  incident  shock  velocity  to  Mach  number  (see  Table  III)  and  reading 
the  calibration  curve,  the  dashed  line  of  Figure  9,  shows  good  agreement  with  the  initially  set 
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Incident  Shock  Temperature 


pressure  ratio  P41  -  as  it  should  for  the  shock  tube  used.  The  departure  between  ideal  and 
actual  P41  at  that  Mach  number  implies  a  smaller  value  of  T5,  in  agreement  with  the 
equilibrium  temperature  of  the  mixture.  This  Martin  Marietta  calibration  curve  was  used  as  a 
standard  of  agreement  and  as  a  basis  of  deducing  the  heat  bath  temperature. 

Ignition  induction  for  diborane  is  presented  on  Figure  16.  Note  that  data  of  Skinner,  et  al 
(Reference  23)  was  obtained  in  a  3-inch  diameter  single-pulse  tube  patterned  after  the  work  of 
Click,  Squire,  and  Hertzberg  (Reference  28),  the  driven  section  being  12  feet  long  and  the  driver 
variable  between  6  and  28  feet.  Thus,  the  L/D  =  48  for  Monsanto’s  driven  tube  is  about  the 
same  (L/D  =  40)  as  Martin  Marietta’s,  indicating  nearly  equal  incident  shock  wave  attenuation 
characteristics.  The  present  measurements  appear  to  extend  the  data  of  Skinner  et  al  by  2  orders 
of  magnitude  with  approximately  the  same  activation  energy. 
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The  activation  energy  obtained  for  points  indicated  on  Table  ill  by  a  method  of  least 
squares  calculation  was  17.5  kcal  compared  to  21.1  kcal  reported  by  the  Monsanto  group  for 
E.R,  =  1.  Note  that  if  Martin  Marietta  were  to  replot  their  data  in  agreement  with  the 
equilibrium  temperature  as  noted  above,  but  not  in  agreement  with  the  calibration  curve,  Martin 
Marietta  activation  energy  would  be  closer  to  Skinners,  for  then  the  40°K,  disagreement  at  the 
high  temperature  end  would  be  removed,  the  disagreement  (from  equilibrium)  at  the  low 
temperature  being  close  to  zero.  However,  if  there  is  a  decrease  in  the  activation  energy  at  high 
temperatures  one  may  conclude  that  it  is  due  to  a  possible  change  in  mechanism  between  the 
high  and  low  temperature  regimes  such  that  the  system  tends  to  be  less  reaction  controlled  at 
the  high  temperature. 

The  efficiency  of  Monsanto’s  tube  is  compared  to  Martin  Marietta’s  and  to  the  ideal  He/Ar 
shock  tube  performance  on  the  expanded  scale  shown  on  Figure  17  (expanded  from  Figure  9). 
The  apparent  better  efficiencies  obtained  by  Monsanto  are  attributed  to  their  use  of  nitrogen 
gas  to  tailor  the  interface  (acoustical  impedance  matching).  However,  Martin  Marietta’s  greater 
departure  from  the  ideal  is  not  a  prima  facie  cause  of  the  slight  departure  from  the  activation 
energy  curve  established  by  the  earlier  investigators. 


Figure  17.  Helium/Argon  Performance  -  Calibration 
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B.  DETONATION  LIMITS 

Detonation  studies,  Figure  18,  give  further  insight  to  differences  between  the  greater 
volume  tube  (Monsanto)  and  Martin  Marietta’s  1-inch  lower-volume  tube.  When  attempts  were 
made  by  the  Monsanto  group  to  measure  ignition  induction  times  of  I^H^/C^  at  temperatures 
greater  than  640° K  for  an  equivalence  ratio  of  1  in  a  99  percent  diluent  mixture,  detonation  of 
the  gas  mixture  occurred  behind  the  incident  shock  wave.  The  driver  pressure  in  these 
experiments  was  4400  mm  Hg.  These  results  showed  that  although  the  gas  mixture  behaves 
dynamically  like  pure  argon  (sec  Figures  14  and  15),  complete  reaction  of  the  1:3  mixture 
theoretically  generates  enough  heat  to  raise  the  gas  temperature  by  about  240°K  at  initial 
pressures  of  about  450  mm  Hg.  This  would  similarly  occur  al  about  45  min  Hg  (Figure  18)  for 
the  present  studies  as  the  detonation  curve  is,  proportionally  speaking,  near  its  peak  there.  Note 
that  the  present  study  is  conducted  at  initial  pressures  that  arc  less  than  the  Monsanto  study  by 
about  one  order  of  magnitude  and  with  a  driver  pressure  of  800  mm  Hgso  that  shock  strengths 
may  be  comparable  and  the  temperature  rise  due  to  combustion  will  be  less  than  240°K, 
proportional  to  the  enthalpies.  The  tailored  interface  technique,  greater  volume, and  high  initial 
pressures  combined  to  exceed  the  explosion  limits  and  to  enhance  the  onset  of  detonation  at 
greater  temper  mrcs,  limits  the  amount  of  ignition  induction  data  obtainable  by  the  large 
volume  tube  at  the  high  temperatures  by  as  much  as  two  orders  of  magnitude  in  the  induction 
time.  On  the  other  hand,  the  low  volume  tube  was  high  pressure  limited  and  data  could  not  be 
obtained  with  sufficient  accuracy  at  the  low  temperatures.  Departure  in  die  calibration  base 
curves  (argon  only)  are  due  to  acoustical  impedance  mismatch  of  the  driver  gases.  To  account 
properly  for  all  differences  it  would  be  necessary  to  extrapolate  data  to  infinite  volume. 


Figure  18.  Comparison  of  Martin  Marietta  and  Monsanto  Deton 
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TABLT  IV 


Wavelengths  Identified.  Spectrum  I 


System 

O 

X  ~  A 
(Publisiicd) 

Measured 

1 

! 

Order 

kef 

Relative  Intensity 

BO-/3 

258 1. 60 

(84.01) 

3 

2 

5.0 

O 

26-1 1 .53 

(40.18) 

2 

i 

5.5 

BO  (3 

2713.80 

( f  5. 1  o ' 

2 

2 

2.0 

^-Kaplan 

2740.00 

(40.60) 

3 

2 

i  .5 

NO  -  X 

2763.30 

(04.00) 

3 

2 

1 .0 

BO  (3 

2999.60 

(97.51; 

3 

2 

>1.5 

OH 

3063.60 

( 64 .  i  1 ) 

3 

2 

5.5 

BH 

3098.90 

(96.10) 

2 

2 

>1.0 

Ar 

3165.30 

!  (63.50) 

2 

i 

1.5 

Ar+ 

3393.7  5 

(94. 10) 

i 

2.0 

BN 

3809.30 

(11.10) 

i 

2 

1.5 

C+ 

3920.68 

(20.78) 

i 

1 

2.0 

B 

4039.57 

(39.75) 

i 

1 

10.0 

N+ 

4133.61 

(33.65  j 

2 

1 

1.0 

co2 

4137.60 

(38.44) 

2 

2 

2.4 

N2 

4141.80 

(4  1 .00) 

2 

2 

2.0 

BO  -a 

4143.40 

(43.03) 

i 

2 

1 .5 

N;< 

4165.00 

(6s. 00)  min 

1,  2 

i 

2.0 

0+ 

4185.456  ; 

(85.63) 

1 

i 

6.0 

10bo? 

4305.89 

(04.88) 

1 

3 

5.0 

BH  “ 

4331.60 

(34 .65)  max  ; 

2 

2 

2.5 

()+ 

4416.97 

(18.10;  1 

i 

2 

7.0 

N+ 

4441.99  i 

(41.10) 

i 

i 

2.0 

1 1  Bl  h 

4502.51 

(01 .10) 

i 

3 

2.5 

CN  “ 

4606.1  0 

(04.78) 

i 

2 

2.5 

Ar 

4872.73 

(7  1.61  ) 

i 

i 

5.5 

Ar 
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(1  7.07) 

1 .  2 

i 

4.0 

,0bo2 

4951 .05 

i 

1 

3 

4.5 

References:  1 .  (I.  K  .  Harrison.  M.l .  I  ,  Wavelength  Tables,  Wiley.  1  956. 

2.  R.  W.  B.  Pearse  and  A.  CJ.  Cavuon,  The  Identificatitm  of  Molecular  Spectra. 
Third  lid..  Chapman  and  Hall.  1963. 


3.  J .  W.  (3.  Johns.  "Tin  Absorption  Spectrum  of  BO->”,  Can  |  of  Physics  39. 
pp  1  738  68,  1  ‘if » I . 


The  third  order  spectrum  on  the  same  photographic  plate  was  also  analysed  and  was 
found  to  yield  the  hydroxyl  (0,0}  bandhead  — n  transition  line  at  3064  A  to  fair 

precision  (1.7  percent  of  30  A)  and  with  55  percent  of  maximum  intensity.  OH  was  not 
observed  in  first  order  because  it  was  obscured  by  the  noise  level  there.  In  second  order  the 
OH  line  was  not  sufficiently  intense  to  be  obtained  below  the  100  percent  transmission  level 
of  the  record.  It  is  to  be  noted  that  OH  emission  has  not  been  observed  in  steady  flames 
during  the  slow  burner  combustion  of  stoichiometric  l^2^6^2  mixtures.  References  23  and  43. 
However,  emission  evidently  exceeded  absorption  by  OH  during  some  or  most  of  the  19 
repeat  runs  producing  the  spectrum  under  analysis. 

A  grey  gas.  Reference  41,  analysis  was  performed,  using  an  absorption  coefficient  of  8.6 
x  10*  cm  “/mole  (Reference  42)  for  OH  at  temperatures  in  excess  of  1100°K.  The 
calculations  show  that  the  absolute  emission  intensity  at  1 160  K.  is  about  0.001  Watts/cm 
while  at  3480  °K  it  increased  to  about  50  Watts/cm2.  We  conclude  that  in  the  detonation 
mode,  the  1<  cal  temperatures  may  reach  sufficiently  large  values  to  promote  large  populations 
of  excited  OH  so  that  emission  exceeds  absorption  as  observed  in  the  present  experiment. 
The  presence  of  OH,  even  in  third  order,  tends  to  qualitatively  verify  both  the  prediction  of 
early  mechanism,  postulated  on  page  19.  and  the  equilibrium  composition  shown  on  Figure 
20. 


Generally,  the  present  observations  are  found  in  agreement  with  emission  spectra 
obtained  from  diborane/air  and  diboranc/NO  flame  studies  by  Wolfhard,  et.  al.,  Reference 
43,  and  oxy  diborane  flame  studies  by  Skinner’s  Group,  Reference  23.  For  example,  the 
intermediate  species  BH,  BO-a  ,  (1  and  B02  were  observed  in  all  the  flame  studies.  Hydroxyl 
radical  was  observed  in  the  noil-stoichiometric  flames.  However,  the  range  of  (non-detonative) 
temperatures  in  which  OH  absorption  dominates  over  emission  by  excited  OH  is  still  to  be 
determined.  This  is  desired  so  that  one  may  minimize  the  guesswork  in  proceeding  with 
quantitative  absorption  measurements.  At  an  equilibrium  temperature  of  1000  °K  in  the 
non-detonative  mode  OH  will  probably  absorb,  according  to  present  calculations  and  to  dat  a 
of  References  44  and  45:  but  if  the  OH  concentrations  arc  as  low  at  1000  "K  as  indicated  on 
Figures  15  and  20,  the  absorption  measurements  may  be  limited  by  instrument  sensitivity. 
Reference  46  It  is  important  to  note  that  the  actual  non-equilibrium  OH  concentrations  to 
be  measured  will  be  one  to  two  orders  of  magnitude  greater  than  those  shown  by  the 
equilibrium  calculations  on  Figure  20.  The  OH  absorption  experiment  was  designed  and  is 
discussed  in  Appendix  A. 
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I.  BACKGROUND 


Research  resulting  from  the  desire  to  utilize  the  processes  of  heat 
release  at  high  power  density  levels  has  greatly  enhanced  the  understand¬ 
ing  o:  phenomena  in  high  speed  gas  dynamics.  Thus,  problems  in  :ombus- 
tion  lynamics  and  reentry  physics  have  received  increasing  attention  by 
vari'.  us  specialists  who  are  finding  a  common  field  of  interest  called  gas- 
wave-d'  namics  (References  47  and  48),  This  new  field  depends  on  experimen¬ 
tal  n  etrods  and  theoretical  developments  in  spectroscopy,  chemical  kinetics, 
statistical  turbulence,  and  radiative  energy  transfer.  Investigative  methods 
are  cemical,  both  in  combustion  dynamics  and  reentry  physics. 

In  the  p  esently  sponsored  AFOSR  investigation  of  high  energy  fuel 
chemistry'  soecial  attention  is  directed  to  the  further  understanding  of 
gas-wave-dynamics  in  supersonic  combustion.  Applications  of  interest 
in  supersonic  combustion  include  development  of  the  supersonic  combus¬ 
tion  ramjet  and  the  use  of  external  burning  for  increased  maneuvering 
capability  of  tactical,  strategic,  and  manned  reentry  vehicles. 

Twc  possible  modes  of  supersonic  combustion  have  been  identified 
or  study,  c.g.,  supersonic  diffusion  flames  and  detonations  (Reference  49). 

The  fo-Tiier  is  largely  deflagrative  burning,  characterized  by  the  mixing 
requirement;  whereas  the  latter  assumes  the  reactive  components  are 
prt  r  xed,  giving  rise  to  the  simplest  explanation,  in  terms  of  Chapman- 
Jo.  _u  t  intonations,  that  the  local  Mach  number  of  the  burned  gas  behind 
the  detoi  ation  wave  is  unity.  A  third  mode  of  supersonic  combustion  is 
pot  sible  that  in  which  the  supersonic  diffusion  flame  achieves  the  con¬ 
dition  of  a  standing  (stationary)  detonation  wave  in  the  reactive  flow  as 
t  is  swe-t  downstream  (References  49,  50). 

Tht  n*  :essity  for  including  the  exact  kinetic  processes  in  the  burning 
region,  in  order  to  understand  the  various  modes  of  supersonic  combus¬ 
tion,  b  is  led  investigators  to  depend  upon  better  known  mechanisms  such 
as  the  branched  H2/O2  reaction  (References  50,  51).  Thus,  a  better  under¬ 
standing  if  the  structure  of  stationary  waves  has  evolved  (References  15  and 
52  thmu  h  56).  It  is  presently  suggested,  for  example,  that  ignition  delay 
in  the  st  itionury  (stable)  detonation  of  H2  is  related  to  the  distance  between 
tne  lu  -n.ng  region  and  the  shock  wave  created  by  pressure  pulses  from  this 
region  .0.,  one  concept  of  a  detonation  wave  which  has  gained  much  ac- 
eeptarv  <  is  that  of  the  shock  wave  followed  by  a  burning  region.  However, 


4.1 


the  study  of  standing  detonation  waves  where  a  free  boundary  is  used, 
such  as  in  the  work  of  Nicholls,  and  detonations  in  tubes  have  shown  some 
interesting  differences.  Whereas  in  a  free  boundary,  the  detonative  flow 
may  be  largely  of  the  "laminar"  type,  in  a  tube  the  reactive  flow  following 
a  running  detonation  wave  becomes  "turbulent"  giving  rise  to  a  non  one¬ 
dimensional  spin  mode  (References  56,  57,  and  58). 

The  origin  of  spin  appears  to  be  related  to  the  scale  of  irregularities 
(turbulence)  in  the  burning  region.  It  is  opined  that  the  turbulence  scale 
depends  on  the  particular  reactants,  the  necessary  activated  state  being 
achieved  in  a  manner  different  (perhaps  more  efficient)  than  in  the  "laminar" 
case.  As  this  turbulence  scale  is  increased,  transverse  acoustical 
waves  coalesce  because  of  superposition  at  the  walls,  giving  rise  tc  the 
shock  waves  observed  within  the  detonative  transition  region  (References 
59  and  60).  Since  these  irregularities  are  not  prima  facie  random,  perhaps 
due  to  anisotropy  induced  by  the  mean  flow,  the  strong  waves  enhance, 
rather  than  cancel  the  motion  creating  an  angular  momentum,  or  span 
(Reference  61).  In  the  case  of  free  boundary  waves,  the  transverse 
acoustical  disturbances  are  present  due  to  combustion  but  are  free  to 
propagate  out  of  the  region  of  influence;  this  may  also  account  for  the 
"laminar"  structure  of  a  standing  detonation  wave. 

There  appears  to  be  some  uncertainty  as  to  the  most  efficient  mode  of 
supersonic  combustion  for  a  given  application.  However,  it  does  seem 
likely  that  attempts  to  utilize  standing  detonation  waves  in  external 
burning  may  require  a  better  understanding  of  the  scale  of  irregularities 
due  to  the  presence  of  the  boundary.  In  any  case,  the  concepts  of  ignition 
induction  and  the  acceptance1  of  shock  tube  kinetics  as  a  means  of  establish¬ 
ing  both  reaction  mechanism  and  ignition  induction  have  been  established. 
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II.  REACTION  PROFILE  MODEL 
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I 


In  the  study  of  reaction  rates  and  mechanisms  of  complex  systems,  it 
is  expedient  to  have  a  means  to  determine  the  reaction  profile,  i.e.,  a 
method  whereby  the  time  rate  of  change  of  each  reacting  species  may  be 
evaluated.  Such  an  analysis  requires  the  simultaneous  solution  of  the  gas 
dynamical  equations  of  change  and  the  species  production  equations.  In  the 
case  of  shock  tube  studies  the  equations  of  change  reduce  to  the  Rankine- 
Huganiot  relations,  and  in  the  case  of  competitive  reactions  a  series  of 
simultaneous  differential  equations  describing  the  relaxation  processes 
must  be  solved.  Hence  the  kinetic  rate  constants  and  the  mechanisms  of 
the  reactions  must  be  known  or  estimated.  In  addition,  a  knowledge  of  the 
coupling  of  the  chemical  kinetics  with  other  relaxation  processes  such  as 
vibration  should,  in  principle,  be  known  since  often  it  is  found  that  the 
chemical  reaction  rates  are  affected  by  changes  in  the  internal  energy 
modes  of  the  species.  Generally  such  effects  are  ignored  initially  until 
there  is  good  reason  to  verify  these  effects  experimentally.  If  it  can  be 
demonstrated  that  the  vibrational  relaxation  time  is  of  the  order  of  or  less 
than  the  chemical  relaxation  time  for  example,  then  provision  must  be 
made  to  account  for  this  effect.  Often  coupling  between  different  relaxation 
modes  may  be  neglected  if  the  temperature  is  not  too  high.  In  addition  to 
these  effects,  temperature  and  pressure  variations  within  the  reaction  zone 
will  also  affect  the  reaction  rates;  however,  the  pressure  variation  is  usu¬ 
ally  small  ai  d,  in  shock  tube  work,  the  temperature  in  the  reaction  zone  is 
very  nearly  constant  if  sufficient  diluent  is  used. 

Duff  (Reference  51)  and  Duff  and  Davidson  (Reference  62)  have  given  a 
method  whereby  the  reaction  profiles  may  be  computed  for  nondetonative 
shock  tube  work  or  the  effects  of  estimated  rate  constants  and  proposed 
reaction  mechanisms  may  be  evaluated.  The  gas  dynamical  equations  and 
the  equation  of  state  are  solved  subject  to  the  initial  conditions,  i.e.,  the 
temperature  and  pressure  behind  the  shock  prior  to  reactions,  the  concen¬ 
tration  of  all  species  at  the  same  point,  and  a  function  of  pressure  and  den¬ 
sity  relabel  to  the  shock  velocity. 


1  h<  .  a  .  :»  s  i  i  !  lii  ratr  <  •  -i >.~>i :  ms  ma\  1  > <  obtained  from  1 1n  i  ito rat  un 
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property  uf  .1  sju  i'M'S  i  as  the  reaction  proceeds.  From  those  data  a  theo- 
rotioal  ear\i-  coniammg  ail instable  para meto rs  may  bo  fitted  to  the  ex¬ 
perimentally  observed  data  and  the  rate  constants  found  from  the  fitted 
values  of  the  parameters  (Heferenee  63).  A  more  sophisticated  method 
is  to  develop  an  algebraic  function  of  the  measured  property  which  is 
linear  with  time.  Values  of  the  property  can  then  be  taken  from  the  data 
at  various  intervals  and  plotted  against  time  so  that  the1  slope  can  be  re¬ 
lated  to  the  rate  constant  (References  64  and  65). 

Duff*  s  reaction  profile  model  was  constructed  for  inviscid  flow,  account¬ 
ing  for  chemical  kinetics. 

Subsequent  formulation  is  aimed  at  the  more  sophisticated  approach; 
following  Marrones'  (Reference  66)  extension  of  Duff' s  idea,  i.e,,  to  in¬ 
clude  vibrational  relaxation  processes  and  chemical  kinetics,  with  their 
cross  couplings,  the  next,  step  is  to  include  the  effects  of  collisional  trans¬ 
port  properties. 
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III.  MATHEMATICAL  FORMULATION 


The  objective  is  to  develop  the  theory  of  planar  shock  wave  propagation 
in  a  reactive  medium,  taking  into  account  the  structure  that  may  arise  due 
to  interaction  effects  of  chemical  kinetics  with  collisional  transport  pro¬ 
perties  and  energy  transfer  mechanisms.  These  effects  superpose  on  the 
general  class  of  inviscid  reaction  profile  models  such  as  Duff' s  and,  more 
recently,  Marrone1  s  (Reference  66).  Planar  detonation  wave  structure 
theories,  which  account  for  the  interaction  effects,  were  investigated  by 
Hirschfelder' s  group  (References  67  and  68),  by  Adamson  (Reference  69) 
and,  more  recently,  by  Petrone  (Reference  70).  In  references  67  and  68 
the  type  of  reactions  studied  was  limited  to  the  exothermic,  irreversible 
and  reversible  unimolecular  reaction  respectively,  with  transport  proper¬ 
ties  grouped  as  constant  values  of  the  Lewis  and  Prandtl  numbers.  Adam¬ 
son  found  an  approximate  analytical  solution  to  the  same  problem,  neglec¬ 
ting  diffusion.  Petrone  considered  a  single  reaction  describing  the  decom¬ 
position  of  a  homogeneous  condensed  explosive  and  included  but  one  trans¬ 
port  property,  namely  artificial  viscosity,  solving  the  time  dependent  one¬ 
dimensional  Lagrangian  equations.  Finally,  a  non-exothermic  reaction 
problem  was  investigated  by  Scala  and  Gordon  (Reference  71)  who  obtained 
an  exact  numerical  solution  to  the  two-dimensional,  time  dependent  Navier- 
Stokes  equations  for  the  supersonic  Row  past  a  blunt  body  where  the  effects 
of  temperature  dependent  transport  properties  and  where  a  system  of  two 
dissociation-recombination  bimolecular  reactions  are  treated. 

In  the  present  development  we  proceed  with  a  general  endothermic  or 
exothermic  scheme,  following  Marrone,  where  the  chemical  source  matrix 
will  be  adapted  to  any  branching  chain  mechanism  and  where  an  explicit 
dependence  of  chemical  reactions  is  introduced  into  the  transport  proper¬ 
ties  subroutines. 

A.  EQUATIONS  OF  MOTION 

The  development  of  chemical  reactions  in  a  branching  chain  mechanism 
may  be  strongly  vibration-dissociation  coupled,  particularly  in  the  initial 
stages  of  thermal  decomposition  during  the  ignition  delay  period.  Thus,  the 
equations  describing  processes  in  the  reaction  zone  downstream  of  a  shock 
wave  must  account  for  such  chemical-mechanical  couplings.  These  relation¬ 
ships  were  discussed  for  a  general  multicomponent  system  in  Reference  72 


and  have  been  reduced  to  a  system  of  first  order,  ordinary,  integro-dif- 
ferential  equations  including  collisional  transport,  radiation  energy  transfer 
and  chemical  source  terms  governing  the  mass  and  energy  transfer  phenomena 
of  a  quasi  steady-state,  one-dimensional,  chemically-relaxing,  multicom¬ 
ponent  gas. 

In  their  simplest  form  these  equations  are: 
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for  energy  transfer.  A,  B,  C,  and  D  arc  constants;  the  e^'s  are  the  reaction 
activation  energies,  a  is  the  Planck  mear  absorption  coefficient,  and  other 
symbols  are  defined  in  Appendix  C. 

Diffusion  mass  flux  and  velocities  Vj  of  chemical  species  are  related 
to  the  concentration,  pressure,  and  temperature  gradients  according  to 


j.  =  p.  v  =  m.  n  v.; 

i  it  tti 


(5) 


that  is,  the  diffusion  velocity  is  given  by  (Reference  72): 
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B.  CHEMICAL  EQUILIBRIA  AND  RELAXATION  PROCESSES 


The  specific  enthalpy  for  chemical  species  is 


h.  *  {'  *  W1)  R,T 


(?) 


where  e-  is  the  specific  energy  content  per  unit  mass  and  the  bracketed 
quantity  is  the  reduced  enthalpy  function 


0<T) 


h(T)  =  y(T) 
pTp  y(T)»  1 


(8) 


a  temperature  dependent  function  of  the  ratio  of  specific  heats,  y,  suggested 
in  Reference  73  for  applications  where  the  state  of  the  gas  may  deviate 
from  thermal  equilibrium.  Considering  Equations  (7)  and  (8)  the  reduced 
enthalpy  is 


/3(T)  =  1  +  £  e  (9) 

f 

where  Cf  is  a  dimensionless  specific  energy  content  per  degree  of  freedom, 
f,  of  the  atomic  or  molecular  system  of  interest.  For  example,  the  energy 
distributions  for  ideal  gases  in  equilibrium  are  described  by 


Thus, 


(10) 


(11) 


where  f  =  3  and  5  respectively  for  atomic  and  diatomic  molecules  in  ac¬ 
cordance  with  equipartition. 

1.  Statistical  Considerations 

These  classical  considerations  from  kinetic  theory  are  not  sufficient 
in  the  case  of  very  strong  shock  waves,  nor  in  the  case  of  detonation  waves, 


49 


and  possibly  not  in  the  cast’  of  somo  deflagrations,  especially  those  due  to 
the  burning  of  high  energy  pyrophoric  fuels.  It  is  therefore  necessary  to 
consider  the  microscopic  model  for  the  distributions  of  chemical  species 
among  the  various  energy  levels  excited  in  the  combustion  processes. 
From  stat  ist  ies, 

•>  •  •'•Tr«r  "2’ 

where  the  partition  function  is  given  by  the  product  of  individual  par¬ 
tition  functions, 

Qf  r°V,'V"U''1  jf  (O) 

associated  with  the  translational,  rotational,  vibrational,  elect  ramie  arid 
nuclear  excitations  of  chemical  species.  Each  partition  function  is  com¬ 
puted  by  an  expression  of  the  form 


'  *  ''  -K../kT 
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J  ‘ 

where  the  E. .  are  the  j  energy  eigenvalues  of  chemical  species  i,  and  the 
g.  are  the  degeneracies  of  these  eigenvalues  in  the  j*  state. 

For  shock  excited  processes  it  is  known  that  the  translational  and 
rotational  degrees  of  freedom  of  atomic  and/or  diatomic  mixtures  equili¬ 
brate  after  about  S  to  r»  and  10  to  .'ton  inlet-particle  collisions  respectively, 
with  increasing  shock  wave  strength.  1  bus,  it  is  safe  to  say  that  equipar- 
t  i  t ion  hoi d s .  i n  gone  r a  1 ,  fo r  t  ran s  1  a 1 1 < w i .  I  h i  s  i  easily  prov  e  because 
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which  yields,  for  all  molecules, 
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in  agreement  with  the  result  from  Equation  (11).  I  he  solution  to 
Schrod i tiger *s  equation  for  tin  rotational  energy  eigenvalues  of  a  rigid 
diatomic  rotator  is 
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and  the  degeneracies  are 
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for  a  rigid  rotator.  Equation  (19)  is  most  easily  summed  by  the  substitution 
of  integration  for  summation,  see  Reference  74, 
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Equations  (12)  and  (13)  then  yield 
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again  in  agreement  with  Equation  (11).  These  results  are  useful  provided 
the  vibrational  levels  of  a  diatomic  gas  are  not  appreciably  excited.  When 
the  vibrational  contribution  is  included, 
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where 

u(T)  =  w-fi/kT. 

With  increasing  shock  strength, 


equation  (22)  reduces  to 
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for  the  "fully"  equilibrated  classical  diatomic  molecule.  Corrections  for 
the  case  of  a  nonrigid,  vibrating  rotator  are  given,  for  example,  in  Refer¬ 
ences  74  and  75.  Equation  (22)  then  becomes 
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where  a  and  ft  are  spectroscopic  parameters  of  the  chemical  species.  Here 
the  similarity  to  classical  kinetic  theory  ends. 


Indeed,  it  may  end  prior  to  the  effects  noticable  when  the  anharmonic 
effect  4>  is  accounted  for.  This  is  because  electronic  states  may  be  excited 
prior  to  the  full  vibrational  equilibrium.  Accounting  for  the  electronic 
excitation,  and  the  heat  of  formation. 
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where  f  is  the  energy  eigenstate  above  the  ground  state  and  where  T  is  the 
equilibrium  temperature. 


2.  Relaxational  Processes 


When  the  state  of  a  molecular  system  is  changed  suddenly,  as  in  the 
case  of  shock  or  detonation  waves,  a  finite  time  is  required  before  the  sys¬ 
tem  reaches  its  new  equilibrium  position.  This  relaxation  time,  r,  is  de¬ 
termined  by  the  efficiency  of  collisional  energy  transfer  between  the  various 
degrees  of  freedom  of  the  system.  The  value  for  t  was  first  obtained  in 
Reference  76.  It  is 
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For  the  present,  consider  only  the  case  of  vibrational  non- equilibrium. 
Then,  the  rate  equation  that  determines  the  deviations  from  equilibrium  is 
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where  TVi  is  the  vibrational  temperature  obtained  by  solving  the  harmonic 
oscillator  relation, 

p  _ _ ^ (30) 

ei  exp^hv/kT^j  - 1 

The  first  term  in  equation  (29)  was  first  given  in  Reference  73  to  describe 
the  rate  of  change  of  vibrational  energy  due  to  collisions.  The  second  term 
accounts  for  the  fact  that  the  average  vibrational  energy  of  the  system  is 
reduced  due  to  dissociation,  while  the  third  term  accounts  for  the  increase 
due  to  recombination.  Thus,  vibrational  and  dissociative  nonequilibrium 
processes  have  been  coupled  according  to  Reference  77. 

Important  vibrational  relaxation  times,  t,  for  high  temperature  air 
species  are 

1  For  Oxygen: 

t  1 
o 


2  For  Nitrogen: 

.  1.11531  xiof  (T,  ,0.5  explm,-]  (32) 

N2  p'  !<t’>  7  J 


t'  ~ seconds  ^ 
pi  ~  dynes/cm 
T'  ~ °K 

C.  EQUATION  OF  STATE 

For  a  real  gas,  the  equation  of  state  is 

p  =  ZpRT,  <33> 

where  Z  is  the  compressibility.  Z  =  1  in  the  case  of  an  ideal  gas  and  will 
be  assumed  so  in  most  instances. 

Equations  (1)  through  (7)  and  (27)  through  (33)  must  be  solved  simul¬ 
taneously  for  the  usual  flow  variables  p(x),  v(x),  p(x),  T(x),  h(x)  and  Xi(x). 


The  species  mole  fractions  X[(x)  obtained  as  a  function  of  the  distance  x 
behind  the  shock  wave,  in  a  one-dimensional  combustion  process,  constitutes 
the  desired  reaction  profile  model. 


D.  KINETICS 

In  the  case  of  combustion  processes  that  are  rate  controlled,  the  dif¬ 
fusion  mass  flux  may  vanish  identically.  Equations  (1)  and  (2)  may  then 
be  rearranged  as  follows: 

N 

dX./dx  -  X.  £m.  dX./dx  <M>  =  N  V /./nv  (34) 

r  ,  i  j  i  r  or 

Equation  (34)  may  then  be  inverted;  as  shown  in  Reference  78,  by  matrix 
techniques  into  a  secular  equation  relating  the  column  matrix  of  mole 
fraction  derivatives,  the  N-  component  snecies  mole  fraction  matrix  of 
rank  N,  and  the  column  matrix  of  N  chemical  source  terms  W.  defined  in 
Reference  72  as  the  mass  rate  of  production  of  species.  1 

E.  SAMPLE  CAT  CUL.  TIONS 

i  \ 

When  oollisional  transport  end  radiative  energy  transfer  mechanisms 
are  suppressed,  the  Martin  Marietta  program  reverts  to  the  CAL  subroutine. 
This  subroutine  rray  be  used  to  compute  the  nonequilibrium  chemical  effects 
due  to  highly'  exothermic,  branching  chain  reactions  during  shock  ignition. 

I 

In  this  case  Equations  (1)  through  (6)  are  altered  by  settings,  k,  D^, 
a  and  equal  to  zero,  differentiating  Equations  (1)  through  (4)  with  re¬ 
spect  to  the  distance  y  following  (he  shock  front,  neglecting  the  radiative 
pressure  term  in  Equation  (3),  and  substituting  for  the  integral  term  in 
Equation  (2)  the  chemical  source  matrix  required.  The  calculation  then 
proceeds  exactly  as  described  in  Reference  Gb. 

The  chemical  system  used  in  this  study  is  the  hydrogen/oxygen  branching 
chain  mechanism.  Table  V  lists  rate  constants  used  for  the  assumed 
mechanism.  These  were  obtained  from  the  work  of  Henrici  and  Bauer, 
Reference  79.  The  result  of  the  seven  species,  eight  reaction  H2/O2 
ignition  induction  calculation  is  shown  on  Figure  23  (solid  line  terminated 
by  circled  points).  Vibrational-dissociation  coupling  was  assumed  un¬ 
important  for  these  calculations.  The  activation  energy  agreed  very  well 
with  that  obtained  by  Hersh,  Frey  and  Gerstein,  Reference  80,  whose  cal¬ 
culations  are  exprapolated  from  the  low  temperature  region.  The  calculated 
activation  energy  seems  to  be  in  slight  disagreement  with  that  measured 
by  Schott  and  Kinsey,  Reference  3.  Departures  in  the  ignition  induction  times 
may  be  explained  in  terms  of  initial  pressure  and  diluent,  Reference  81. 


54 


Ignition  Induction  Time  ~msec 


ami. i  v 


Hair  Constant:;  fun  Assumes)  II, /O^  Mechanism 

Krai’  lion 

Rat t-  Constant : 

( cc  mule'  *  sec”  ' ) 

II,,  t  M  »2  II  t  m 

2.23  x  1012  x  T  1  exp( - 92  600/RT) 

HH 

O,  »  M  -  20  t  M 

1  ft  1 

3.6x10  x  T  exp<-  1  18  OOO/RT) 

pips 

II.,  •  O.,  *  2  011 

2.51  x  !012  exp(  -  39  OOO/HT) 

1 

il  <  O.,  —  OH  t  O 

7.75  x  1013  exp( -  1 4  450/RT) 

O  i  11,,  -  Oil  t  II 

1013  exp(-  8  800/RT) 

0 

OH  *■  11.,  -  II  , O  t  11 

2.3  x  1013  exp( - 5  150/RT) 

f 

Oil  +  OH  -  11^0  +  O 
!!.,0  +  M  -  1!  +  OH  +  M 

1.55  x  1012 

g.h 

5  x  1014  exp(-  105  OOO/RT) 

i 

The  rate  constants  shown  have  been  determined  for  M  is  argon. 

a.  A.  I..  Myerson  and  W.  S.  Watt,  J.  Chem.  Phys.  49,  425  (1968). 

b.  K.  L.  Wray,  in  Progress  in  Astronautics  and  Rocketry,  F.  R.  Riddel, 

Ed.  (Academic  Press  Inc.,  New  York,  1962),  Vol.  7.  p.  181. 

c.  Value  suggested  by  G.  L.  Schott. 

d.  D.  Gutman  and  G.  L.  Schott,  .1.  Chem.  Phys.  46,  4576  (1967). 

e.  I.  M.  Campbell  and  B.  A.  Thrush,  Trans.  Faraday  Soc.  64,  1265  (  1968). 

f.  W.  E.  Wilson,  Western  States  Section  Meeting,  The  Combustion  Institute, 
LaJolla,  Calif.,  1967,  Report  on  the  Establishment  of  Chemical  Kinetics 
Tables,  Chemical  Propulsion  Information  Agency,  April  1967. 

g.  W.  E.  Wilson  and  J.  T.  O' Donovan,  J.  Chem.  Phys.  47,  5455  (1967). 

h.  The  reference  quoted  gives  no  temperature  dependence.  In  earlier  work 
(F.  Kaufman  and  F.  P.  Del  Greco,  Symp.  Combust.  9th  Cornell  Univ., 
Ithaca.  N.  Y. ,  1963,  659  ( 1953))  it  was  suggested  that  the  activation 
energy  cannot  be  more  than  (1-2)  kcal/mole.  In  our  calculations  we  used 
the  room-temperature  value  without  an  activation  energy. 


Calculations  (points)  are  also  shown  for  eight  species,  ten  reaction  and 
nine  species,  eleven  reaction  models  assumed.  These  contained  the  ter- 
molecular  reaction  H  +  Og  +  Ar  -*•  HO2  +  Ar,  and  the  collisional  dissociation 
reaction. 


OH  +  Ar  -*  O+H+Ar. 

Suppressing  the  termolecular  reaction  had  the  effect  of  increasing  the 
ignition  induction  time  as  shown  at  a  nominal  temperature  of  1820°K.  The 
same  calculation,  repeated  at  Cornell  University  produced  a  still  higher 
value  of  the  ignition  delay  time  at  1855JK.  This  is  attributed  to  machine 
differences  since  the  same  model  was  used.  Otherwise,  differences  in 
the  model  used,  though  deemed  adequate  for  ignition  processes  does  tend 
to  make  a  great  difference  in  the  ignition  time  and  may  produce  errors 
in  the  activation  energy.  This  was  discovered  after  performing  a  calcula¬ 
tion  with  the  nine  species,  eleven  reaction  model  of  Reference  82. 
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The  output  lists  species  concentrations  variation  with  distance  down¬ 
stream  of  the  shock  front,  or  as  a  function  of  time  from  the  front.  The 
reverse  of  each  reaction  is  accounted  for  in  the  CAL  subroutine.  For 
these  calculations  the  induction  time  is  defined  as  the  time  for  maximum 
OH  concentration.  The  (final)  equilibrium  temperature  occurring  in  the 
reaction  is  used  to  obtain  the  ignition  induction  plot.  Figure  24  shows  the 
OH  concentration  envelope  for  all  reactions  from  which  the  maximum  OH 
concentration  was  obtained  at  2.35  /j sec,  for  an  equilibrium  temperature 
of  1855°K.  Figure  25  gives  the  positive  Q.  .  (i  =  OH)  matrix  elements  which 
represent  the  instantaneous  net.  rates  of  formation  of  species  i  in  reaction 
j.  The  envelope  of  the  absolute  values  of  all  such  matrix  elements  for  OH, 
viz  Figure  24,  was  used  to  obtain  the  ignition  induction  time. 
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Figure  25. 


Rate  of  Production  of  OH,  QoH.j  vers^s  Time  of  Several 
Reactions  in  the  ilydrogen/Oxygen  Chain 
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IV.  EXTENDED  MODEL 

Of  immediate  interest  to  the  Martin  Marietta  Corporation  is  the  chemi¬ 
cally  coupled  shock  wave,  reaction  profile  model  now  coded.  This  model 
is  constructed  with  arbitrary  endothermic  and  exothermic  chemical  source 
input.  It  includes  as  principal  subroutine,  the  quasi-steady  state,  one¬ 
dimensional,  inviscid,  reacting  air  chemistry  model  of  Reference  66.  As 
such,  this  subroutine  is  ready  for  use  in  high  temperature,  chemical  non¬ 
equilibrium  problems,  both  endothermic  and  exothermic,  provided  that 
radiative  energy  and  collisional  transport,  effects  are  unimportant.  A 
number  of  practical  engineering  problems  presently  require  the  CAL 
(subroutine)  inviscid  model  to  supplement: 

1  Propulsion, 

2  Heat  shield,  and 

j?  Communications 

systems  design  of  advanced,  high  acceleration  vehicles. 

The  CAL  subroutine  is  used  as  a  diagnostic  device  for  determination 
of  reaction  rate  data  within  the  framework  of  self  consistent,  branched- 
chain  reaction  mechanisms  (Reference  79).  To  better  achieve  such  eval¬ 
uations,  Reference  66  is  being  extended  to  include  the  interacting  effects 
of  finite  rate  chemical  reactions,  radiative  energy,  and  collisional  trans¬ 
port  properties  of  high  temperature  multicomponent  systems.  To  the  de¬ 
gree  that  collisional  transport  models  are  known,  Reference  83  for  example, 
for  purposes  of  calculating  the  multicomponent  viscosity,  diffusion,  and 
thermal  conductivity  coefficients,  the  proposed  extension  is  complete.  How¬ 
ever,  even  the  simplest  radiative  energy  transfer  model  proposed,  in  the 
optically  thin  limit,  will  be  of  limited  use  due  to  the  paucity  of  photon  ab¬ 
sorption  coefficient  data,  particularly  in  the  chemiluminescent  regions  of 
interest.  Recently  there  has  been  a  great  demand,  and  corresponding  ac¬ 
tivity,  to  obtain  these  data,  along  with  refinements  of  the  radiative  energy 
transfer  models,  both  at  the  optically  thin  and  thick  limits  and  within  the 
framework  of  a  precise  molecular  band  structured  model,  see  Reference  84. 

It  is  believed  that  the  extended  model  will  play  a  significant  role  in  de¬ 
sign  problems  of  interest  to  the  Nation's  aerospace  engineering  effort  be¬ 
cause,  in  addition  to  the  determination  of  chemical  nonequilibrium  effects, 
it  will  also  be  possible  to  estimate: 


i 
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1  Radiative  heating  due  to  the  rocket  motor  and  the  heat  shield 

2  Effects  of  viscosity,  diffusion,  and  thermal  conductivity,  for  applica- 
tions  where  the  nondimensional  transport  property  parameters 
(Schmidt  and  Prandtl  numbers)  are  variable 

3  The  generation  and  populations  of  free  photoelectrons  that  contribute 
to  plasma  noise  effects,  disrupting  communication  signals. 


A.  EQUATIONS  CODED 

It  is  desired  to  investigate  the  interacting  effects  of  multicompc  nent 
collisional  transport  properties,  radiative  energy  transfer  and  chemical 
kinetics.  It  is  known.  Reference  83  for  example,  that  collisional  momentum 
and  mass  transfer  are  locally  independent  of  chemical  kinetics;  but  depend 
on  the  previous  history  of  physico-chemical  processes  in  the  flow  through 
the  local  multicomponent  molecular  weight  and  temperature.  In  a  reacting 
flow,  particularly  where  supersonic  combustion  may  give  rise  to  steep  tem¬ 
perature  and  concentration  gradients,  it  is  important  to  access  this  mutual 
interaction.  Collisional  energy  transfer  is  far  more  complicated  than  the 
other  two  transport  properties  because  of  internal  energy  exchanges  be¬ 
tween  species  occupying  varied  excited  energy  states.  Reference  85  to  92, 
in  addition  to  the  energy  transport  dependence  on  chemical  kinetics  and 
prior  history  of  processes.  These  combined  energy  transport  phenomena 
are  grouped  into  the  coefficient  of  thermal  conductivity. 


Starting  with  general  equations  applicable  to  chemical  reactions  in  flow 
systems  we  obtained  a  system  of  first  order,  ordinary  integro-differential 
equations  accounting  for  collisional  transport,  radiation  energy  transfer,  and 
chemical  source  terms  governing  transfer  phenomena  in  a  quasi-steady 
state,  one-dimensional,  chemically  (and  vibrationally)  relaxing,  multicom¬ 
ponent  gas.  Thermodynamic  expressions  accounting  for  the  contributions 
of  excited  vibrational  and  electronic  states,  including  anharmonicity,  have 
been  incorporated  into  the  energy  equation,  whereas  Landau-Teller-Hertzfeld- 
Treanor  and  Marrone  theory  is  used  to  account  for  coupled  vibrational- 
dissociative  nonequilibrium,  (refer  to  Equation  1-33),  see  References  66 
and  93. 

In  their  present  form,  the  equations  have  been  further  developed  to  adopt 
Treanor's  modification  to  the  l;unge-Kutta  method  (Reference  94).  A  matrix 
of  Runge-Kutta  coefficients  has  been  programmed  from  the  following  set  of 
17  equations,  and  machine  logic  has  been  developed  for  the  overall  program: 
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(44) 
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^  =  U  (50) 

dy 

dU  =  d2U  =  _3_dp  +  3pU  ^  (51) 

dy  dy2  "4a  dy  4p 


These  equations  are  used  to  compute,  in  order  of  their  appearance,  the 
rate  of  change  (with  respect  to  distance  y  from  the  shock  front)  of 


r. 

j 


(MW)  = 


X  L— 

the  vibrational  energy  of  jin  species,  calories/mole; 

y  -u./u,  concentration  of  jtl1  species,  moles/unit  mass,  including 
the  ^effect  of  multicomponent  diffusion  through  the  ratio  of  species 
velocity  Uj  to  mass  averaged  velocity,  u; 

mole  fraction  of  jth  species; 

total  molecular  weight  at  a  point  y  downstream  of  the  shock  front 
due  to  the  production  and  recombination  of  species,  including  the 
effect  of  the  rate  of  change  of  their  concentrations,  gms/mole; 
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Y  .  -  concentrations  of  species,  moles/unit  mass; 

j.  =  mass  flux,  gm/cm^s; 

p  =  mass  density,  gm/cc; 

p  =  pressure,  dynes/cm^; 

T  -  temperature,  degrees  Kelvin; 

T  =  rate  of  change  of  temperature  with  respect  to  y; 

U  =  rate  of  change  of  mass  averaged  velocity  with  respect  to  y, 

B.  VIBRATION  -  DISSOCIATION  COUPLING  MODEL 

The  nonpreferential  vibrational  energy  model,  Equation  (35),  is  derived 
in  Reference  93  and  is  included  in  the  principal  (C.A.L.)  subroutine,  where 

X.  =  the  vibrational  relaxation  distance  of  j**1  species,  cm; 

Oj  =  the  characteristic  vibrational  temperature  of  j**1  species,  °K; 

I'm  s  Ty"1  -  T^”1  is  a  temperature  parameter,  °K  \  relating  the  vi- 

J  brational  and  translational  temperatures; 

i  u 

Nj  =  the  number  of  vibrational  energy  levels  of  j  species; 

iL  t  h 

A^  =  1  or  zero,  is  the  vibrational  coupling  of  j  species  to  the  i 
J  reaction  depending  upon  which  reaction  will  be  affected  by  the 
vibration-dissociation  coupling  process; 

Q. .  =  is  the  mole-volumetric  rate  of  production  of  jth  species  from  the 
ith  reaction; 

X  =  is  the  degree  of  nonequilibrium  of  the  itl1  reaction. 

In  the  nonpreferential  model,  dissociation  is  hypothesized  to  occur  with 
equal  probability  from  any  level  of  vibrational  excitation  of  jtfl  species. 
However,  recent  theoretical  arguments,  References  95  and  96,  suggest 
that  dissociation  is  coupled  preferentially  to  the  upper  vibrational  levels 
of  excited  molecules.  Thus,  a  preferential  vibrational  energy  model  is 
also  included  in  the  CAL  subroutine.  Presently,  such  models  exist  only 
for  diatomic  molecules. 

C.  COLLISIONAL  TRANSPORT  PROPERTIES 
1.  Mass  Transfer  -  The  Diffusion  Coefficients 

Equation  (37)  requires  an  expression  for  calculating  the  multicomponent 
diffusion  coefficient,  D  .  For  this  prupose,  Reference  83  gives 
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where  ||  K  ||  is  t he  determinant  of 


K. 
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(MW) 

_  J 
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r  (MW). 

1  1 

(52) 


(53) 


and  L  are  dimensional  velocity  and  scale  parameters  in  the  CAL  sub¬ 
routine,  the  K0^  are  the  minors  of  ||k||,  the  K..  =0,  and  the  first  approxi¬ 
mation  to  the  binary  diffusion  coefficient  is  computed  from  the  formula 


*  (0.002  628) 


r 

(MW).  +  (MW) 

J  0 

1  2< 

:MW).  (MW) 
_ 1 _ 9L _ 

2(1,1)  * 
12  °12  (T12  ) 


(54) 


Equation  (54)  shows  that  the  binary  diffusion  coefficient  varies  directly 
with  the  3/2  power  of  the  local  translational  temperature  and  varies  in¬ 
versely  with  the  local  pressure,  the  square  of  the  collision  diameter  a12 
and  the  collision  integral  1  (T^*).  where  T12*  *  kT/c^  is  the  reduced 
temperature,  e^/k  is  the  potential  parameter  in  °K,  and  k  is  Boltzmann's 
constant.  Equation  (54)  was  checked  out  for  self  diffusion.  Table  VI. 


TABLE  VI 


Comparison  of  Calculated  and  Observed 
Self  Diffusion  Coefficients  at  1  atm* 


Gas 

T(  °K) 

2 

cm 

Present 

-1 

5 

Ref 

83 

Experimental 
(cm^  sec"1) 

AR 

353.2 

wsm 

0.249  +  0.003 

273.2 

IffiK 

0.156  +  0.002 

0.158  +  0.002 

77.7 

0.0136 

0.133 

0.0134"+  0.0002 

N„ 

353.2 

0.273 

0.273 

0.287  +  0.009 

273,2 

0.174 

0.174 

0.185  +  0.006 

0.172  +  0.002 

77.7 

0.0161 

0.0161 

0.0168-  +  0.0003 

*Check  1  P581  using  Equation  (3.2-46)  Reference  83. 
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Collision  diameters  and  integrals  arc  tabulated  in  References  83  and 
07  for  like  particle  encounters  (single  species  gas).  The  collision  diame¬ 
ters  are  averaged  for  pairs  of  interacting  (unlike)  molecules,  according  to 
Reference  83, 

=  ^■<ct1+<72)  (55) 


where  the  a  -  are  obtained  from  Tables  I- A,  page  1110,  Reference  5  and 
Table  1,  Reference  07.  The  required  collision  integrals  are  obtainable 
as  a  function  of  T^*.  where  the  potential  parameter  for  pairs  of  interact¬ 
ing  molecules  is,  again  according  to  Reference  83, 

f12  =  (r>6) 

2.  Momentum  Transfer  -  The  Viscosity  Coefficients 

For  a  single  species  gas  the  first  approximation  to  the  viscosity  co¬ 
efficient,  from  Reference  83,  is 


/T(MW). 

M,  x  10°  =  26.693— (57) 

a.  J2U»^(T*) 


is  tabulated  extensively  for  molecular  species  of  interest  in  Reference  21. 
The  first  approximation  to  the  multicomponent  viscosity  coefficient,  re¬ 
quired  in  Equations  (49)  and  (51),  is  given  by  Reference  83; 
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where  the  diagonal  elements  are 
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the  off-diagonal  elements  are, 
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the  collision  integral  ratio  is 
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(61) 


tabulated  in  Tables  I-N  and  VII-E  (pages  1128  and  1176  respectively)  of 
Reference  83  for  the  Lennard- Jones  and  Buckingham  potentials;  and  the 
(molar)  universal  gas  constant  is 


R 


8.313404  x  107 


ergs/°K/mole, 


3.  Energy  Transfer  -  The  Thermal  Conductivity  Coefficients 


Reference  83,  page  1196,  recommends  use  of  the  formula 

K  =  l"K  •  1  +  +  K  ( 62 ) 

L  rnixjj  int  react 

with  which  to  estimate  thermal  conductivity  due  to  a  multicomponent  mix¬ 
ture  of  polyatomic,  polar  and  nonpolar  molecules  in  chemical  nonequili¬ 
brium,  and  under  flow  conditions  where  local  thermodynamic  equilibrium 
(l.t.e.)  holds;  i.e.,  where  collisional  and  radiative  processes  do  not  appre¬ 
ciably  alter  the  distribution  of  molecular  velocities  from  the  Maxwell- 
Boltzmann  distribution  at  the  local  bath  temperature. 


The  first  term  on  the  right  hand  side  of  Equation  (62)  is  the  first  approxi¬ 
mation  to  the  coefficient  of  thermal  conductivity,  not  including  the  effects  of 
internal  energy  transfer  between  excited  states  occupied  by  interacting 
molecules,  nor  the  energy  transfer  effects  of  chemical  kinetics,  but  account¬ 
ing  for  the  transformation  of  internal  energy  into  translational  motion;  i.e., 
for  a  multicomponent  gas. 
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where,  in  the  nomenclature  of  Reference  66,  the  diagonal  elements  are. 
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the  off-diagonal  elements  are. 
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and  the  first  approximation  thermal  conductivity  coefficient  for  a  single 
component  gas,  including  the  effect  of  internal  energy  exchange  with  the 
translational  modes,  is 

[-j],  =i[I5-6« 
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-  (15- 


/  (MW). 
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(66) 


The  first  term  on  the  right  hand  side  of  Equation  (66)  represents  the 
first  approximation  to  the  thermal  conductivity  for  polyatomic  gases  where 
it  was  assumed  that  the  internal  energy  of  a  molecule  does  not  depend  on 
tne  molecular  velocity,  and  that  the  exchange  of  internal  and  translational 


energy  occurs  sufficiently  rapid  so  that  l.t.e.  holds  at  each  point  of  the  flow. 
The  second  term  represents  the  correction,  for  nonpolar  gases,  by  Mason 
and  Monchick  (Reference  (86) 


where 

(cv)k  is  the  specific  heat  of  the  kth  internal  mode  and 

Zk  is  the  number  of  collisions  required  to  exchange  a  quantum 
of  energy  of  this  mode  with  the  translational  degrees  of 
freedom. 


Experimental  data  for  a  number  of  common  gases  were  examined  in 
Reference  86  where  it  is  concluded  that  better  agreement  between  experi¬ 
ment  and  theory  occurs  at  temperatures  less  than  about  400°K  when  Equa¬ 
tion  (66)  is  used  to  correlate  data.  However,  at  high  temperatures  two 
compensating  effects  take  place  in  that  the  translational  conductivity  is 
less  than  the  value  predicted  by  the  first  term  on  the  right  hand  side  of 
Equation  (66)  while  the  reciprocal  Schmidt  number. 
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is  ijivnU’i'  than  unity;  i.e.,  greater  than  the  Eucken  correction,  Reference  83, 
which  was  used  prior  to  t ho  theoretical  prediction  by  Mason  and  Monchick. 
Further,  since  only  a  few  of  the  large  number  of  relaxation  times  needed  to 
evaluate  Z|.  are  known,  a  reciprocal  Schmidt  number  value  of  £  =  1.32  is 
used  in  the  present  work  in  order  to  use  the  data  computed  for  the  first  term 


of  Equation  (32)  tabulated  in  Reference  97.  The  actual  form  of  Equation  66 
used  for  monatomic  gases  check  out  is  Equation  8.2-31  of  Reference  83, 
see  Table  MIL  The  Eucken  correction  was  used  for  polyatomic  gases, 
Fable  VIII.  Equation  63-66  were  used  to  check  out  the  binary  gas  mixtures 
shown  on  Figure  2  6  and  Table  IX,  An  option  to  compute  the  exact  value  of 
[Kj]l  *s  included  in  the  program  for  those  cases  where  the  relaxation  times 
have  been  measured  (diatomic  gases)  and  to  be  consistent  with  the  relaxa¬ 
tion  phenomena  calculation  included  within  the  framework  of  the  Vibration- 
Dissociation  Coupling  Model  (the  CAE  subroutine)  above. 


TABLE  VII 
Thermal  Conductivity 

Thermal  Conductivity  of  Monatomic  Gases:  Comparison  of  Experimental  Values 
with  those  Calculated  Using  Equation  (8.2-31)* 


7  1-1 

X  10  in  cal  cm  sec  deg 


Helium 

— 

Neon 

Ar 

gon 

Kryston 

Xenon 

mm 

Calcu- 

Ex  peri- 

Ggmm 

Expert- 

Calcu- 

gjgjmgi 

Ex  peri- 

mm 

1  at  ec 

mental 

mental 

lated 

mental 

90.2 

496 

489 

140 

141 

68 

1 

40 

III 

194.7 

870 

876 

297 

293 

149 

152 

85 

WBm 

273.2 

1097 

1097 

398 

385 

206 

190 

120 

3406 

1110 

390 

208 

124 

3438 

394 

3510 

373.2 

4329 

4165 

1340 

1357 

509 

506 

271 

272 

161 

168 

•<91.2 

5173 

4947 

1589 

1595 

622 

614 

339 

340 

204 

208 

379.1 

5752 

5504 

1770 

1789 

699 

685 

384 

388 

233 

237 

•Check  l  P573  Reference  83 
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TAlfl.E  VIII 


I'lu- final  Conduct ivity 


Thermal  Conductivity  of  l*o|yatom  10  Oases:  Comparison  of  Experimental 
Values  with  those  Calculated  Using  Formula  (8.2-33)* 


\ 

107  in  c 

il  c:m  St 

-  1  . 

•c  (leg 

1 

r 

Hydrogen 

Oxyg 

on 

Carbon 

Dioxide 

Methane 

Nitrous 

Oxide 

°K 

Calcu- 

Ex  peri- 

Calcu- 

Ex  peri- 

Calcu- 

Expert- 

Calcu- 

Ex  pen- 

Cal  eu- 

Kxjx'ri- 

luted 

mental 

lated 

mental 

lated 

mental 

lated 

m ental 

lated 

menial 

100 

1580 

1625 

224 

216 

2  54 

254 

150 

327 

321 

200 

2094 

3064 

435 

438 

183 

227 

508 

522 

456 

425 

2  73 

3871 

3965 

568 

584 

267 

349 

696 

734 

578 

567 

4040 

577 

360 

720 

300 

4099 

4227 

614 

635 

300 

398 

764 

819 

618 

619 

*Clu’el 

2  P574  Reference 

83. 

table:  ix 

Thermal  Conductivity 

Comparison  of  Calculated  and  Experimental  Values  for  Binary  Gas  Mixture  Using  Equation 
(8.2-42)  and  Mix  Formula  on  Page  1196: 


X  10 

5  .  ,  -1  -1  ,  -1 

m  cal  cm  sec  deg 

Helium  -  Argon 

Hydrogen  -  Carbon  Dioxide 

Hydi  ogen  -  Argon 

Percent  of 

1  ightpr 
Constituent 

Experimental 

Calculated 

Percent  of 
Lighter 
Constituent 

Experimental 

Calculated 

Percent  of 
Lighter 
Constituent 

Experimental 

Calculated 

0.0 

3.80 

3.08 

0.0 

3.60 

3.91 

0.0 

3.90 

3.99 

27.04 

7.42 

7.43 

10.0 

5.10 

5.10 

9.0 

5.50 

5.58 

4  5.37 

10.77 

10.81 

14.2 

5.70 

5.67 

18.0 

7.30 

7.37 

8  '  .08 

23.20 

24.43 

25.0 

7.70 

7.33 

40.0 

12.60 

12.88 

04.61 

20.30 

30.84 

35.5 

10.00 

9.30 

60.0 

18.70 

20.00 

100.00 

33.86 

35.24 

50.0 

13.50 

12.85 

80.2 

27.00 

30.87 

75.0 

22.70 

23.10 

100.0 

40.00 

4H.99 

90.0 

31.50 

34.97 

100.0 

_ 

40.40 

L 

48.49 

♦Check  3  P576-578  Reference  83. 
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Thermal  Cm 


Figure  26.  Check  Calculation  of  Thermal  Conductivity  for  Various  Binary 

Systems 
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The  second  term  on  the  right  hand  of  Equation  (62)  is,  according  to  Ref¬ 
erence  83,  page  1196  and  Reference  87, 
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(68) 


c  is  the  specific  heat  at  constant  pressure  of  the  jth  gas  component 
Pj  (species)  and 

. .  is  the  coefficient  of  self  diffusion  of  the  jth  species,  derivable 
^  from  Equation  (54). 


According  to  Butler  and  Brokaw,  Reference  90,  who  have  shown  that  for 
a  set  of  s  chemical  species  [m]  participating  in  r  chemical  reactions  which 
can  be  written  stoicniometrically 


s 

JH  v..  [Mj  =  0;  i  =  1.2 . r  (69) 

j=l  1J 


it  follows  that  the  third  term  of  Equation  (62)  is 
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where 


(71) 


is  the  enthalpy  change  of  the  ith  reaction  and  the  Hj  are  the  species 
enthalpies  referred  to  a  common  base.  The  determinant  elements  are 
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(72) 
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If  the  chemical  reactions  are  slow,  the  thermal  conductivity,  Equation  (70) 
depends  on  the  boundary  conditions,  see  Reference  83  (page  1197),  Refer¬ 
ences  91  and  92. 


D.  RADIATIVE  ENERGY  TRANSFER 

The  scope  of  present  work  includes  investigation  of  gas  dynamically 
coupled  radiative  energy  transfer,  with  particular  emphasis  on  applications 
in  supersonic  combustion.  Such  coupled  interactions  are  included  within 
the  framework  of  the  reaction  profile  model,  the  sources  being  due  to 
chemiluminescence  and  inelastic  collisional  encounters  among  the  reactants. 


Beginning  with  the  simplest  approximation,  that  of  the  optically  thin  gas, 
the  energy  conservation  term  expressing  the  divergence  of  radiative  flux 
in  Equation  (49)  is  (see  References  98  and  99) 
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where 
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(74) 


is  the  Planck  mean  absorption  coefficient. 
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.hv/kT^ 

is  the  Planck  function  which,  integrated  over  all  frequencies  gives 


(75) 
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The  Stefan- Roltzman  constant  is 


(77) 


(78) 


and  the  spectral  intensity  is 

I  =  k  B  .  (79) 

V  V  V 

To  evaluate  the  Planck  mean-absorption  coefficient,  combine  Equations 
(  74)  -  (  78)  to  obtain 
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5.6697  x  10"  erg/cm  °K  s; 


the  absorption  coefficient  is 
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v  v 
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(l-e-h^T). 


t3e'4d(A(t) 

7T  Sk 


(80) 


where 


C  =  hv/kT  (81) 

and 

*(U  =  (pr,  v)”1  (82) 

is  the  radiative  mean-free-path. 

Two  alternate  methods  of  determining  absorption  coefficients  (e.g.,  re¬ 
ciprocal  Equation  (82))  or  spectral  intensities  (Equation  (79))  are  available 
to  evaluate  the  required  integral,  Equation  (88).  The  first  method  requires 
a  direct  experimental  measurement  of  1/*  over  the  wavelength  range  of 
interest,  at  the  equilibrium  temperature.  Such  measurements  are  not 
available  for  the  special  case  of  interest;  i.e.,  the  combustion  of  pyrophoric 
fuels.  The  second  method  depends  on  quantum  mechanics.  In  this  case, 
measurements  of  oscillator  strengths,  transition  frequencies  and  spectral 
constants  are  still,  for  the  most  part,  required.  Further,  as  is  pointed  out 
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by  Herzberg,  Reference  100,  a  knowledge  of  the  excited  states  of  interme¬ 
diate  species  is  of  great  importance  in  the  determination  of  chemilumines¬ 
cence  and  in  understanding  the  occurrence  of  elementary  processes;  i.e., 
the  elementary  chemical  kinetic  steps  in  the  reaction  mechanism.  These 
problems  may  be  approached  as  follows. 

1.  Local  Thermodyanie  Equilibrium 

Chandrasekhar's  hypothesis  of  local  thermodynamic  equilibrium,  Ref¬ 
erence  101,  is  an  important  concept  applied  in  the  calculation  of  gas  dynam¬ 
ically  coupled  spectral  intensities  and  radiative  energy  transfer.  Accord¬ 
ing  to  this  hypothesis,  at  each  point  (p,  T)  of  the  flow  field,  one  defines  a 
gas  at  temperature  T  and  specie  density  p  such  that  p  is  the  population  of 
excited  energy  states  of  that  gas  in  true  thermodynamic  equilibrium  at 
temperature  T.  That  is,  the  state  populations  are  locally  defined  by  a 
Boltzmann  distribution,  exp  (-K/kT).  The  condition  of  local  thermody¬ 
namic  equilibrium  does  not  require  that  the  radiation  field  be  in  equilibrium 
with  the  flow  field  nor  is  chemical  equilibrium  a  requirement  for  thermody¬ 
namic  equilibrium  on  a  local  basis. 

A  subroutine  code  for  calculation  of  the  spectral  intensities  of  species 
in  nonequilibrium  high  temperature  gas  is  required.  The  reaction  energies 
and  populations  of  continuum  states  obtained  from  the  CAL  subroutine 
must  he  taken  into  account.  The  following  section  gives  a  review  of  the 
basic  formulation  and  recently  accepted  approximations  for  calculation  of 
spectral  intensities  of  diatomic  species  in  either  equilibrium  or  chemical 
nonequilibrium  in  the  optically  thin  limit. 

2.  Review  of  Basic  Formulism 

The  radiation  process  in  a  model  sample  of  high  temperature  gas  may 
be  regarded  as  the  emission  of  photons  of  spectral  frequencies  v  <:/ X 
caused  by  electronic  (orbital)  transitions  of  a  molecule  from  an  initial  state 
of  energy  K  to  a  final  state  of  energy  Em:  where  En  >  Km  is  the  condition 
that  radiation  of  frequency  (Em  -  Kn)/h  -  vnm  be  emitted.  The  rate  of 
emission  of  photons  of  energy  hvtl(ll  from  a  single  source  in  the  sample  is 
called  the  radiation  intensity  If  (watts/molecule).  Classically,  the  inten¬ 
sity  of  the  radiation  field  is  determined  by  Fourier  analyzing  the  spectrum 
due  to  a  dipole  source.  For  a  single  (orbital)  electron,  the  time  averaged 
inte  nsity  corresponding  to  the  #»h  Fourier  component  of  frequency  tv 
(accounting  for  two  polarizations  of  the  field)  is 

\t  --  ^  {2* tv  )4  !pj  |2/ c'*  (83) 
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where  is  the  Fourier  coefficient  of  the  dipole  moment.  Now,  of  all  mole¬ 
cules  in  the  sample  which  occupy  the  initial  state,  let  Anm  be  the  fraction 
which  participates  in  the  transition  from  state  n  to  state  m  per  second. 

Then,  on  a  statistical  basis 


I.  =  A  hv 
i  nm  nm 


(84) 


Comparing  Equations  (83)  and  (84)  obtains  the  Einstein  coefficient;  i.e.,  the 
transition  probability 
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nm 
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64tt  /  nm! 
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(85) 


where  Fourier  terms  of  Equation  (83)  have  been  formally  replaced  by  the 
matrix  elements  of  the  dipole  moment;  e.g.. 


<m  |p 


k|n>  =  ^  ^  ^*(m)  pk  4<n)  dr 


(86) 


The  state  eigenfunctions  i|j  are  obtained  by  solving  the  time-independent 
Schroedinger  equation,  given  the  appropriate  potential  function.  One  must 
sum  Equation  (85)  over  all  states  m  which  have  energy  Em  <  En  to  obtain 
the  total  probability  per  unit  of  time  that  state  n  will  undergo  spontaneous 
depopulation  through  emission  of  light  photons.  The  fact  that  spontaneous 
transitions  incur  isotropic  radiation,  while  induced  emission  (such  as  in 
Raman  spectra)  is  aligned  with  the  incident  field,  leads  to  the  Einstein 
probability  Bnm  that  the  system  will  undergo  a  induced  emission  from  n 
to  m  under  the  influence  of  an  external  field.  On  the  other  hand,  induced 
absorption  m-*n  is  also  likely;  thus,  accounting  for  possible  degeneracies 
gk  in  the  electronic  states  of  a  molecule,  the  Einstein  coefficient  for  in¬ 
duced  absorption  is, References  102  and  103. 
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mn 


(87) 


Since  the  intensities  of  molecular  electronic  transitions  from  an  upper 
state  n  to  a  lower  state  m  are  determined  by  the  transition  probability 
Anm  in  emission  or  Bnm  in  absorption,  one  can  theoretically  find  the 
emission  intensity  from  Equations  (84),  (85)  and  (86)  after  multiplying  (84) 
by  the  number  density  Nm  of  molecules  populating  the  lower  electronic 
state.  However,  the  absolute  intensity  of  electronic  transitions  is  usually 
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determined  experimentally  from  the  absorption  spectrum  because  of  the 
general  difficulty  of  finding  Nm.  The  development  chosen  here  is  the  ab¬ 
sorption  spectrum  approach  to  facilitate  inclusion  into  the  program  of 
new  measurements  of  physical  data  that  appear  in  the  literature. 

Now,  the  spectral  intensity  due  to  an  optically  thin  slab  of  gas  is  given 
approximately  by  Kirchoff' s  law;  viz, 

\  -  BA  (88) 

where  k  is  the  absorption  coefficient,  and 

Bx =  ~  ^”ir  [exp<hc/kTK>  -  l]~l  <89) 

4  ir  \ 

is  the  Planck  energy  density  rate  per  steradian.  In  that  hc/kT  is  large, 
Wien's  approximation  is  good.  Therefore,  Equations  (88)  and  (89)  com¬ 
bine  to  give 


2hc2k. 


>X  * 


exp(-hc/kT\). 


(90) 


Given  the  absorption  coefficient  k\  relating  the  energy  density  ID  (incident 
per  unit  time  on  a  "one- dimensional"  sample  of  gas)  to  the  transmitted 
energy. 


-k  Ax 

I  =  I  e  X  (91) 
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one  integrates  over  the  pertinent  bandwidth  interval  to  obtain  the  absorp¬ 
tion  intensity.  Reference  104, 
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Therefore,  the  integrated  absorption  coefficient  is 
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where  Nm  is  the  number  density  of  molecules  populating  the  lower  elec¬ 
tronic  state,  rn  is  the  classical  electron  "radius",  e^/mc^,  and  the  "f- 
number"  is  defined  by  Equations  ( 85),  (87),  and  (93).  In  terms  of  the  matrix 
elements  of  the  dipole  moment 

2 

mv 

f  _  _ nm 

T  '  2 
3hc 

Equation  (60)  is  the  "electronic  oscillator  strength"  f,  a  dimensionless 
quantity  that  can  be  found  experimentally  and  is  most  important  to  the  cal¬ 
culation  of  both  transition  probabilities  and  spectral  intensities.* 

Differentiating  Equation  (59)  yields  the  absorption  coefficient 

k.  *  N  ffr  c  (95) 

X.  m  o  dv 

Penner,  Reference  106,  approximates  the  derivative  in  the  following  man¬ 
ner.  Assuming  all  transitions  from  lower  energy  levels  are  to  the  lowest  r 

upper  energy  level,  the  derivative  is  equal  to  the  total  oscillator  strength 
times  the  rate  of  change  of  the  fraction  of  molecules  that  are  in  an  excited 
energy  state 


<m|pk  |n> 


(94) 


Thus, 


i.e 


•  < 


(96) 


*Bethe  and  Salpeter,  Reference  105,  shows  that  for  hydrogen-like  atoms  of 
nuclear  charge  Z  «  137,  the  dipole  approximation  holds  and  2^  f  ■  1. 
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The  present  formulation  uses  an  approach,  originally  applied  to  NO  absorp¬ 
tion  spectra.  Reference  107  which  is  in  essential  agreement  with  Penner. 
That  is,  employing  the  rotating  oscillator  model  for  diatomic  molecules, 
the  rate  of  change  of  the  excited  population  varies  directly  with  the  Franck- 
Condon  electronic  transition  probability  for  a  single  vibrational  mode  with 
average  rotational  overtones.  Thus,  introducing  rotational  and  vibrational 
partition  functions  Qr",  Qv"  in  the  lower  (double  primed)  state.  Keck,  et  al.. 
Reference  108,  show 


ill 

dv 
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Substituting  Equation  (97)  into  (95)  yields  the  absorption  coefficient  for  a 
particular  vibrational  transition  and  averaged  rotational  structure.  Thus, 


Jr  r  f  N  e 
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oo 


Qr"  Qvi, 


B*  -  B"  I 
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-(E  ,  +  E  .  )/kT 

Vv',e  r  V  (98) 


A  summation  of  Equation  (98)  over  all  vibrational  transitions  which  can 
occur  at  the  wavelength  X  will  give  the  absorption  coefficient  at  X  for  the 
electronic  transition  involved: 


-hc(v_  -  v)/kT 

-<Er.  +  Eyl)/kT 
(99) 

A  close  inspection  of  the  exponential  term  Er»  will  reveal  the  simple  ex¬ 
pression  for  the  manner  in  which  the  summation  is  to  be  executed!.  Con¬ 
sider  a  vibration  band  whose  rotational  structure  degrades  to  the  red,  i.e., 
greater  X.  The  summation,  in  this  case,  should  be  over  all  vibrational 
bands  for  which  the  band  head  Xv«v"  is  equal  to  or  smaller  than  the  wave¬ 
length  in  question  (Xv«v"  s  X),  Reversing  the  case  for  vibration  bands 
which  degrade  to  the  violet,  the  summation  applies  to  ail  band  heads  equal 
to  or  greater  than  the  wavelength  in  question  (Xviy"  a  X).  Further,  it  may 


#Most  symbols  not  defined  are  the  accepted  notation:  e.g.,  v  and  r  for 
vibrational  and  rotational  quantum  numbers, 

tThe  summation  accounts  for  branching  of  the  spectra,  see  Reference  109. 
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be  noted  from  the  spectrographic  constants  involved,  that  the  quantity 
(Bg  -  B”)  will  be  positive  for  all  bands  degrading  to  the  violet,  and  nega¬ 
tive  for  those  degrading  to  the  red.  Substituting  these  facts  into  the  expres¬ 
sion  for  Er«  it  is  easily  verified  that  the  summation  for  both  types  of  de¬ 
grading  is  that  over  all  transitions  for  which  Er'  2  0,  i.e.. 
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(100) 


Finally,  Equation  (99)  may  be  substituted  into  Equation  (90)  obtaining  the 
spectral  intensity 
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A  factor  <$>  is  incorporated  in  the  derivation  of  Reference  108,  which  may  be 
extracted  from  the  above  expression.  In  terms  of  $,  Equation  (101)  is  written 
in  the  compact  form, 

I.  *  2jt  r  he2  4>  X."5  N  f{£exp(-hcv  /  T)  (102) 

X  o  m  kT  oo 


where 
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(103)* 


Rotational  and  vibrational  partition  functions  for  the  lower  state  population 
are 


kT 


he  B' 


(104) 


to  a  good  approximation,  and 


The  notation  Er'  2  0  under  the  summation  sign  in  the  expression  for  $ 
implies  that  all  terms  will  be  added  into  the  sum  except  those  for  which 
Eri  <  0. 
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exactly. 

The  (Fortrat)  averaged  rotational  energy  level  of  the  emitting  state  is 
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where 


v  ,  n  =  v  +  (E  ,  -  E  n)  T  ; 
v'v  oo  v'  v  he 


(107) 


the  vibrational  energy  of  the  emitting  state  is 


E  ,  =  he  ( v  |  “  v  ,). 
v 1  v 1  o 1 


(108) 


Assuming  the  rotating  oscillator  model  for  diatomic  molecules. 


vv.  s  <v'  +  1/2)  -  x'e  (v*  +  1/2)2  +  y l  (v*  +  1/2)3  (109) 
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Likewise,  for  the  absorbing  or  lower  state. 


(110) 
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V  =  W"  (v"  +  1/2)  -  u/J  x"  (v"  +  1/2)2  +  w"  y”  <v"  +  1/2 )3  (112) 


(113) 


Including  the  Free-Bound  Radiation,  e.g.,  Referencea  110  and  111 
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where  a  is  the  photodetachment  cross-section,  cm  . 


Note  that  spectral  intensities  calculated  with  this  subroutine  can  be  used 
to  predict  the  optically  thin  limit  radiation  signature  of  shock  layers  and 
trails  that  are  predominantly  composed  of  O2,  and  NO  in  chemical 
equilibrium  (or  in  chemical  non  equilibrium  when  the  populations  of  excited 
states  are  known). 
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APPENDIX  A 


OPTICAL  ABSOKBIMKTRY  OF  THE  B2II6/02  SYSTEM 

Specific  radiation  absorbimetry  measurements  are  to  be  made,  with  mi¬ 
crosecond  resolution,  to  monitor  the  appearance  and  disappearance  of  the 
hydroxyl  *OH)  radical,  postulated  as  a  cnain  carrier  intermediate  in  the 
high  temperature  oxidation  of  diborane  (B2Hg). 

The  reactor  to  be  used  is  a  one-inch  diameter  chemical  shock  tube  with 
helium  gas  driving  shock  waves  in  argon,  generating  temperatures  from 
500  to  1000°K  in  the  region  behind  the  reflected  shock  wave.  The  source 
of  OH  radiation  to  be  used  is  a  flash  lamp  containing  water  vapor.  The 
source  radiation  absorption  will  be  measured  in  the  homogeneous  gas  sam¬ 
ple  behind  the  reflected  shock  wave  in  argon  at  initial  pressure  near  76  mm 
Hg  and  containing  the  reactants,  B2Hg  and  02,  in  amounts  of  the  order  of  1 
percent  by  volume. 

These  measurements  will  give  the  change  in  the  concentration  of  OH  with 
time.  From  this  it  may  be  possible  to  deduce  the  effective  activation  ener¬ 
gy  and  corresponding  reaction  rate  (s)  leading  to  a  better  understanding  of 
the  high  temperature  reaction  kinetics  of  B2H^  oxidation.  It  may  also  be 
possible  to  verify  the  postulated  chain  mechanism  following  ignition  induc¬ 
tion. 

Measurement  of  the  OH  concentration  is  to  be  made  by  ultraviolet  light 
absorption,  using  a  flash  lamp  as  the  source  of  characteristic  OH  radiation. 
The  OH  radical  was  selected  for  this  study  because  of  its  known  relative, 
high  stability  and  very  high  optical  absorption  coefficient  (about  2,3  x  10® 
cm*/ mole).  The  integrated  absorption  coefficient  for  OH,  over-  a  selected 
number  of  rotational  lines,  will  be  measured  by  comparing  the  observed 
absorption  with  the  computed  OH  concentration  at  equilibrium. 

Using  an  analytical  technique  developed  at  Cornell  University  -  beginning 
with  Bauer,  Schott  and  Duff  (1958),  the  concentration  of  a  specific  chemical 
species  as  a  function  of  time  at  high  temperatures  may  be  obtained  for 
species  with  fine  line  absorption  spectra,  such  as  the  OH  radical.  This  can 
be  done  by  using  as  a  source  an  intense  fine  line  emission  spectrum  gener¬ 
ated  by  exiting  the  species  w'hich  one  wishes  to  observe  in  the  sample. 
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Ill  this  case  the  source  of  Oil  radiation  is  a  pulsed  discharge  through  a 
short  capillary  containing  water  vapor  at  a  pressure  of  about  10  mm  Hg 
(Figure  27).  The  lamp  is  viewed  through  the  ring-shaped  anode.  The  elec¬ 
trical  circuit  for  this  flash  lamp  consists  of  a  3.7  microfad  capacitor 
charged  to  10  kilovolts  in  series  with  a  1.5  kilohm  resistor  and  an  8  milli¬ 
henry  choke  (Figure  28).  This  cirauit  produces  a  light  pulse  which  rises 
in  about  300  microseconds  and  's  nearly  constant  for  the  following  1  milli¬ 
second,  during  which  absorption  measurements  are  made.  The  entire  ex¬ 
perimental  setup  is  shown  schematically  on  Figure  29. 


C4H10  10  mml  _ 

N2  10  mm/ 

CS2  20  mm 

Ar  10  mm/ 

Figure  27.  Characteristic  Radiation  Source  Lamp 
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30  mh  Each 


Output:  5  kV  Across 
Secondary  to  Start  Breakdown 


Figure  28.  Schematic  of  High  Voltage  Capacitor  Bank  and  Lovotron  Switch 
to  be  Used  for  Energizing  the  Characteristic  Radiation  Source 


Lovotron 

Switch 


Figure  29.  Absorptivity  Measurement  Instrumentation 
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APPENDIX  B 


A  SHOCK  WAVE  ANALYSIS 


One-dimensional  conservation  equations  were  used  for  computation  of 
flow  variables  in  the  chemical  equilibrium  region  downstream  of  shock 
waves  traversing  a  multicomponent  gas  mixture. 


A.  ONE -DIMENSIONAL  CONSERVATION  EQUATIONS 


These  equations  are: 


P1U1 


P2U2 


for  miss  conservation. 


P1  +  P1U1  ’  p2+  P2U2 


for  momentum  conservation,  and 

hi  +iui2  *  h2  +iu22  •  (3) 

for  energy  conservation.  The  terms  p,  u,  p,  and  h  are  the  gas  density, 
velocity,  pressure  and  specific  enthalpy  designated  upstream  (subscript  1) 
and  downstream  (subscript  2)  of  the  shock  wave.  Since  the  specific  enthalpy 

h  *  e  +  p/p  (4) 

introduces  a  fifth  variable,  namely  the  internal  energy,  e,  per  unit  mass 
of  gas  component,  essentially  a  function  of  temperature,  the  system  of 
equations  to  be  solved  requires  equation  (4)  in  order  to  calculate  enthalpy 
and  an  equation  of  state, 

p  =  pRT  =  NkT,  (5) 

to  calculate  temperature.  N  is  the  total  number  density  of  v- component 

species,  N  *  2  N-  and  the  density  p  =  2  m,N{. 

i-1  i=l  1 

Rearranging  equation  (3)  and  introducing  equation  (1) 

1  2  ..  .  .  1  ,  ,  ,2  2 

2  “i  "  *h2  “  hl*  +  2  Pl/P2  U1 
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then 


hence 


|uj2  [l  -  <p/p2>2]  *(h2-hi> 

Uj2  [l  -  (p/p/  ]  ■  2<h2  -  h2)  . 


(6) 


Rearranging  equation  (2),  introducing  equation  (1) 


2  2 


thus 


then 


plul  ’(p2'pl)  +  p2(pl/p2)  U1 


u,2  *  <P  -  Pl)/P,  +  (Pl/P2) 


“l2  [*  -  pl/p2]  =  (P2  ‘  pl)/pl  ' 

Dividing  equation  (6)  by  equation  (7) 

1  +  p/p2  =  2pl(h2  "  hl)//{p2  "  Pl) 

then 

<h2  -  hl) 


(7) 


i  +  p,/p2  -  ^o/p,)  ,p2/Pi  .  ,) 


(8) 


Introducing  the  equation  of  state  and  its  ratio  from  region  1  to  region  2 
across  the  incident  shock  wave  one  obtains  from  equation  (5)  and  equation 
(8) 

2(h  -  hJ/R.T 
i  +  p,/p2  •  2  *  * -1- 

r  /  \ a 2  i 

(p2'pVBjT1  "  1 

2(h2  -  h,) 


(9) 


R2T2*p2,/pl)  '  R1T1 

From  equation  (5)  the  gas  constant  is  defined  (for  equilibrium  temperature 

T)  as  v  v 

r/ S  m  n.  -  k^.Z  N.  (10) 

^ i= 1  1  l  i«l  i 
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i.e.. 


R  =  N  k  £  N./  2  N.M.  =  N  k/  2  X.M. 


°  i=l  1  i=l  1  1  °  i=l 


l  l 


then 


R  =  N  k/M  =  8.3143  x  10?/M  ergs/°K/gm  (11) 


where  the  molecular  weight 


m  =  2  X.M. 

1=1  1  1 


(12) 


is  obtained  by  averaging  the  v-component  molecular  weights  over  their 
mole  fraction  X.  =  N./N.  (13) 


i  i 

Equation  (9)  becomes 


i  +  p/p2  - 


2<h2  *  hl>  M2 


N  k 
o 


T2(p2//pl)  "  T1 


5 

M 


1  - 


(14) 


thus,  solving  for  Pg/Pj 


N0k<T2/ra2» 


T1  M2 

p,/p. - -  —  (1  +  p./p,)  *  2(h  -  h  ) 

T2  “I 


(15) 


where  the  enthalpy  may  be  obtained  from  JANAF  tabulations  according  to 

h2  "  hi  1  (H2/M2  ’  Hi/KV  JE7  (16> 


and  at  T  *  298°K,  1^(298)  a  0.  Therefore,  equations  (15  and  (16)  become 


N  kT 

O  i 


II 

-  T2  Mu 


(1  +  Pl/p2)  *  2H2JE7  .  (17) 


Let  N  kT„  •  or,  then,  expanding  (17): 

O  a 
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1 


M  1  M 

»  *  »  'A  '  NokTl<i^  (p2/p,)  -  NokT.  <17 ' 2H2JE7 


a(p2/px)  + 

therefore. 


M 


a  -  2H_JE7  -  N  kT.^) 
2  o  1  M1 


JVL 

p2/pi  -  NokTi(s:)  *  ° 


P2^P1  = 


5* 

*  -  2H2JE7  -  NCT^) 


/2a 


2  ] 

,Mg.  ' 

“  -  2H2JE7  -  NkT,^ 

i. 

+  4*NokTi(M7j 

where 


M, 

*  ■  «  -  2H2JE7  -  NokTl(M? 


and 


M. 


V  ■  N  kT.d?)  . 
o  1  M, 


(18) 

(19) 

(20) 


Note: 


J  *  4.1840  joules/calorie,  NQk  ■  8.3143E7  ergs/TC/mole,  (21) 


and 


H 


■  Ex/\^+H  ^  calories/mole 


2  .  ,  i  i 

i»l 


(22) 


where  the  h^  are  enthalpies  for  the  individual  V2  stable  components  in 
the  equilibrium  region  downstream  of  the  incident  shock  wave  and  Ho'  ' 
represents  the  bond  energy  or  heat  of  dissociation  of  the  system.  From 
equations  (6)  and  (16),  the  incident  shock  velocity,  is  obtained,  namely. 


“i  y 


r - 

r  1 

/2H2JE7/W2 

' '  W'f  . 

(23) 
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B.  COMPUTATION  PROCEDURE 


The  following  procedures  will  be  used: 

Initial  conditions  are  specified  as  pj  and  T.  =  Tpjyj  (usually  298°); 
vi  components  at  Tj  are  given,  along  with  their  mole  fractions  (or 
partial  pressures).  V2  components  are  specified  according  to 
stoichiometery. 

2  T2  is  the  equilibrium  temperature  required,  with  V2  equilibrium 
concentrations  at  T2  determined  from  a  free-energy  minimization 
subroutine,  following  that  of  Reference  112  and  113. 


a  In  the  computation  scheme,  obtain  the  molecular  weights: 

M  »  2  and  M  =  £  X.(2)M.*2)  , 

1  .  .  1  1  2  .  ,  1  i 

i=l  i=l 


(2) 

where  the  v2  mole  fractions,  X^  ,  are  determined  through 
iteration  required  to  satisfy  the  free  energy  minimization 
(equilibrium  constants)  and  mass  balance  constraints.  The  con¬ 
centrations  may  then  be  computed  according  to 

c.(2)  ■  (x/2^  p  /ffl )E3,  moles/liter  .  (24) 

1  i  Ct  & 


b  Obtain  H2  from  equation  (22),  with 

\<2)*H,-H288 

taken  from  JANAF  thermochemical  tables  or  computed  using  the 
method  of  References  112  and  113  available  as  an  optional  alternate 
subroutine  to  Martin  Marietta  program  shock  (F4>  018). 


£  The  gas  dynamical  quantity  P2  is  obtained  from  equations  (18), 

(19)  and  (20)  and  pi  is  determined  from  the  initial  conditions 
specified.  A  second  iteration  loop  is  required  between  the  main 
program  (Shock)  and  the  equilibrium  subroutine  so  that  p2  is 
simultaneously  satisfied  by  equations  (18)  and  the  v2  equilibrium 
concentrations,  i.e.,  the  final  pressure  p2  obtained  at  the  specified 
equilibrium  temperature  T2  must  satisfy  the  gas  dynamical  rela¬ 
tion  from  equation  (7):  2 


1  - 


o/oy 


(25) 


9! 


to  be  calculated  after  finding  p_/p  and  II  from  equations  (18)  and  (23) 

^  1  1 

C.  OTHER  CALCULATIONS  REQUIRED 
Consider  the  shock  tube  just  before  and  after  diaphragm  burst. 


Before 


11  JJ.L 


A.fter 


E  j — C  I 


(4)  (3)  (2)  (1) 


The  initially  set  pressures  are  pi  for  the  gas  mixture  in  region  (1)  and 
p4  in  the  high  pressure  helium  driver  gas  region  (4).  After  diaphragm 
burst,  the  incident  shock  wave  I  is  driven  to  the  reflection  end  plate  by 
contact  surface  C,  while  an  expansion  wave  E  moves  in  the  opposite 
direction.  Gas  particles  flow  from  (1)  to  (2)  so  that  the  excess  gas  velo¬ 
city  in  region  (3)  is,  according  to  Ref  : 


— 4  V4  ”  1 

U3SU1*U2S^1  1  -  <P3/p4>17T 


ui(1  •  W 


where  equation  (1)  was  used  in  the  last  step. 

Since  the  pressures  are  equal  across  the  contact  surface,  P2  =  P3» 
equations  (26)  and  (27)  may  be  rearranged  to  give  the  diaphragm  burst 
pressure  ratio.  0 


P41  '  P21 


y4  -  1  Uj(l  -  Pj/p2) 


where  equation  (25)  may  be  written 


Pai  -  *  *>»»■. 


A 


The  parameters 


c 


pi 


2  X.(1)e  (1)(T  ) 

.  ,  L  p.  1 

1  =  1  *1 


1  (0 
.2,  X.(c  -  R  ) 
i=l  i  p.  o 

i 


(30) 


and 


M,  '  u/f  1RoT/Ml 


(31) 


are  the  ratio  of  specific  heats  and  the  Mach  number  respectively.  Equa¬ 
tion  (29)  may  be  used  alternately  for  frozen  flow  calculations  (no  chemical 
reactions)  using  the  formulas  given  for  p2/pi»  for  fixed  -y  and  M,  Refer¬ 
ence  114.  The  frozen  flow  temperature,  T2,  will  be  greater  than  the  corre¬ 
sponding  value  of  the  equilibrium  temperature  at  the  same  value  of  the 
incident  shock  velocity. 

Program  Shock,  which  in  part  was  obtained  through  Cornell  University, 
computes  equations  (18),  (23),  (25)  and  (28)  simultaneously  to  determine 
the  diaphragm  burst  pressure  ratio  as  a  function  of  ui  or  Mp  For  ideal 
(frozen  flow)  diaphragm  burst,  the  equation  used  (Reference  )  is 


P41  * P2 / 


1  - 


v4  ~  1 

+  1  » 


',4T4Ml  <M‘ 


M 


1  J 


(32) 


The  equilibrium  calculations  for  flow  across  the  reflected  shock  wave 
follow  an  iterative  procedure  similar  to  the  one  outlined  above.  The  value 
of  (T5)gQ  is  obtained  through  a  second  double  iteration  in  program  Shock 
which  requires  a  stoichiometric  mass  balance  and  free  energy  minimiza¬ 
tion  constraint.  However,  the  pure  helium/argon  gas  couple  used  in  the 
shock  tube  design  requires  no  consideration  of  chemistry  so  that  the  ideal 
shock  relations  were  used  for  chat  purpose. 
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APPENDIX  C 
LIST  OF  SYMBOLS 

denotes  vibrational  coupling  of  species  to  i^h  reaction 
Planck  functions 


velocity  of  light,  cm/s 


X  L 

specific  heat,  at  constant  pressure  of  j  n  species,  Cp  =  c'  /R'* 
(nondimensional)  j  ° 

specific  heat  at  constant  volume  of  j**1  species 

multicomponent  diffusion  coefficient 

binary  diffusion  coefficient 

self  diffusion  coefficient 

electronic  charge 

rotational  energy  of  the  emitting  state  relative  to  the  ground 
(lowest  upper)  state  of  the  electronic  band,  ergs 

vibrational  energy  of  the  emitting  state  relative  to  the  ground 
state  of  the  electronic  band,  ergs 

vibrational  energy  of  the  lower  state,  ergs 

electronic  oscillator  strength,  or  f-number,  dimensionless 

Planck's  constant,  Ti  =  h/2ir,  erg  sec 

enthalpy  of  species  including  heat  of  formation;  (dimension- 

less)  hj  >  hVRJTJ, 

intensity  in  watts/cm'‘/steradian/micron 

species  diffusion  flux  =  p.U. 

J  J 

Boltzmann's  constant 

Planck  mean  absorption  coefficient 

reference  length 


^Primes  denote  dimensional  quantities. 
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M.  a  [M] 

(MW) 

(MW). 

3 

m. 


J 

m  =  pU 

m 

N. 

N 
N 


m 


n. 

J 


Qr„(T) 

Qv„(Tx) 

Q.. 

1J 

qv'v" 


R' 

o 


T 

T 

U 


m. 


V. 

J 

X. 

3 

y 

zk 

r. 

j 


denotes  species 

molecular  weight  of  gas  mixture  defined  by  Equation  (4) 

t  h 

molecular  weight  of  j  n  species 

number  of  electronic  levels  included  for  jth  species 
mass  flow  rate 
electron  mass 

number  of  vibrational  energy  levels  of  j  species 

Avogadro's  number  =  6.023  x  10^  per  AMU 

number  of  molecules  (per  cm^)  populating  the  lower  electronic 

state 

number  of  atoms  per  species 

pressure;  p  =  P'/P^UM  (nondimensional) 

rotational  partition  function  of  the  lower  state 

vibrational  partition  function  of  the  lower  state 

mole -volumetric  rate  of  production  of  M^  from  reaction  i 

Franck-Condon  factor  for  the  transition  from  v'  to  v"  (square 
of  the  overlap  integral),  dimensionless 

universal  gas  constant  (R1  =  1.98647  cal/mole  °K  in  CAL  sub¬ 
routine) 

equilibrium  gas  temperature,  °K 

parameter  of  jth  species  having  units  of  temperature 
velocity;  U  =  U'/U^  (nondimensional) 
vibrational  coupling  factor  of  jth  species 
mole  fraction  of  jth  species 

coordinate  normal  to  shock;  y  =  y‘/L  (nondimensional) 
number  of  collisions  per  mode 

concentration  of  jth  species,  moles/unit  mass  including  effect 
of  diffusion  flux,  Equation  (8) 

concentration  of  species  in  moles  per  unit  mass  (see  Equa¬ 
tion  (6)) 

ratio  of  specific  heats  for  gas 

vibrational  energy  of  jth  species;  ^  (nondimensional) 
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r 

r 

; 

defined  by  Equation  (47) 

r* 

I 

i 

®vj 

characteristic  vibrational  temperature  of  j**1  species 

K 

coefficient  of  thermal  conductivity 

f” 

A 

nondimens ionalizing  term,  see  Reference  4  (the  CAL  subroutine) 

A  =  (U')2MW'/R,T» 

o  o  o  o 

r 

X 

j 

X 

vibrational  relaxation  distance  of  j**1  species 

*  t  ■ 

wavelength 

:  • 

! 

u 

coefficient  of  viscosity 

V 

V 

wavenumber  of  the  vth  vibrational  level,  cm”* 

r 

Vv1' 

wavenumber  for  the  vibrational  transition  between  the  states  v> 

and  v",  cm"* 

V 

wavenumber  of  the  (0,0)  vibrational  transition,  cm”* 

■  ■  ' 

00 

V  .  , 

stoichiometric  coefficients  of  j**1  species  in  i**1  reaction 

1J 

s 

reciprocal  Schmidt  number 

f,’ 

p 

density;  p  =  p'/p*  (nondimensional) 
o 

r 

cr 

o 

photodetachment  cross-section,  cm  j  collisional  diameter 

!,  *  - 

V 

degree  of  nonequilibrium  of  i**1  reaction 

r  r 

to. 

3 

Be  ’ 

vibrational  frequency  of  jtn  molecule 

|| 

•  •  i  : 

'  1  i  ' 

OJ 

e 

>  Spectral  constants.  The  1  refers  to  the  excited  electronic  or 

:  1 

u  X 
p  p 

emitting  state  and  the  M  refers  to  the  lower  electronic  or  ab¬ 

■ 

l 

sorbing  state. 

1 

...  L 

1 

!  '  * 

Subscripts 

I 

I  ' 

00 

refers  to  vibrational  equilibrium 

I 

| 

o 

free  stream  conditions  ahead  of  shock  wave 

1 

i 

pertaining  to  i**1  reaction 

I 

I  ! 

j 

pertaining  to  species 

I 

1 . 

n 

pertaining  to  nt!l  electronic  level 

l 

V 

pertaining  to  vtl1  vibrational  level 

1 

1 

r 

r 

pertaining  to  r^  rotational  level 

1 

i 
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